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EINTRODGCTION 


Every object in our environment emits thermal radiation. The peak energy of this 
radiation 1s characteristic of the object’s temperature and. for temperatures commonly 
found in our environment. is concentrated mostly in the infrared region of the 
electromagnetic spectrum. This physical phenomenon has found many applications 
since the second world war as it enables passive detection of targets by their self- 
emission alone. | 

Passive search, detection. identification and tracking have alwavs been of great 
importance in naval warfare. Also, the necessity to improve our combative capabilities 
in areas where radar performance is degraded (low flving small cross section targets) 1s 
imminent. Infrared technology promises much in this area. The rapid growth of 
technology in solid state integrated circuits, microprocessors, infrared detectors and the 
better understanding of infrared physics during the recent decades has made it possible 
for sophisticated infrared systems to be designed and developed. 

A very important step in designing an infrared optical system is the knowledge of 
the spatial and spectral radiances of the target and its background as well as the 
atmospheric effects on the emitted thermal radiation. Almost always, targets appear 
against a background that complicates the detection process since radiation from 
Objects surrounding the target 1s confused with radiation emitted or reflected from the 
target. In the analysis of thermal imaging systems the target-to-background radiance 
contrast is frequently described by the equivalent tempcrature difference. In the case of 
observation at low elevation angles above the sea surface the surface radiance 1s 
Sraomel iimitenced by the reflection of the radiation emitted by the sky. The actual 
thermometric temperature difference is no longer a good representation of the radiance 
contrast. An “effective temperature difference” can however be defined taking into 
account the reflected radiation. Anv new information and additional data about 
background scenes will improve the understanding of these effects and may be used to 
optimize computer algorithms or to unplement new or improved techniques in signal 
processing to minimize background effects, thus improving system capabilities and 


optimizing performance. 


This thesis describes the collection and analvsis of thermal images of an 
instrumented ship target and its background under various well defined environmental 
conditions, followed by the development of a model for the effective target-to- 
background temperature difference. The target was the oceanographic ship R V 
“POINT SUR” and the background was sea suriace andeskg, Wyhe data weremeollcered: 
stored. and processed using the “AGA 780 THERMOVISION”, a hquid nitrogen 
cooled thermographic device with digital recording and data analvsis hardware. All 
measurenients were taken from a distance of about 650 meters using the Mercury 
Cadnuuni Telturide (IlgCdTe) detector of the Thermovision 780 scanning system in the 
8 - 14 jtm spectral region. Five experiments were conducted during an operational 
oceanography student crurse on nud May 1986 and an adequate number of thermal! 
image pictures was taken for various environmental conditions on certain dates. 

Analysing these thermal images using the measured ship surface emissivitv and 
the AGA software package. both the radiance and the effective temperature 
distributions of the ship were obtained. The atmospheric transmuttance for this purpose 
was coniputed using the LOWTRAN 6 code. and the emissivity was measured from 
samples of the ship paint. The effective sea surface temperature was also obtained from 
these pictures. It was observed that the apparent sea surface temperature was affected 
by reflected skv radiation, and a model was developed to describe the effective target- 
to-background temperature difference under varying environmental conditions. 

Later, three more experinients were conducted in nud Noveinberw$ts6 in ordemame 
obtain further, more accurate, data under different environmental conditions, after the 
Thermovision had been calibrated. This was done to check the previous results and 
iniprove the information on the reflected sky thermal radiation and the effective target- 
to-background temiperature difference relation. 

The present work consists of eisht chapters. [n the first Chapter sweustaccate 
background problem. In the second chapter we discuss the background radiances which 
affected our measurements, the marine background and the sky background. The 
atmospheric propagation problem and the LOWTRAN code used to find the 
atmospheric transmittance for the various environmental conditions are discussed in 
chapter three. Chapter four provides a familiarization with the equipment used to 
collect the data. In the fifth chapter we describe the calibration problem and discuss 
the derivation of the empirical calibration relations. In the sixth chapter we discuss 


how the thermal data we had collected were analvzed to give the temperature and 
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radiance distribution of the target. A further analvsis of these data, to see how the sea 
background radiance was affected bv reflected sky thermal radiation, is given in 
Chapter seven. In this chapter we also calculate the apparent sea surface temperature 
based on environmental parameters and compare the target-to-background temperature 
aiepemce sensed by On detector with the actual thermometric temperature difference 
measurements. In the last chapter we discuss our conclusions and recomendations for 
further work. The discussion and analvsis of the data obtained during the mid 
November measurements are given in Appendix A. In Appendix B the TI-59 program 
used to get the in-band radiant flux for the radiance calculations is explained. 

The results of this work will lead to a better understanding of the background 
effects and may add some information to be used for the implementation of techniques 


in signal processing. 


Il. TARGET AND BACKGROUND 


A. BACKGROUND RADIATION NOISE 

Background radiation noise is nolse generated in an optical detector by photon 
flux reaching it from all points other than the target. This kind of noise usually 1s 
separated into two components [Ref. 1]: 

1. Temporal, which ts due to fluctuations in the background emitted photons. 
2, Spatial, which is due to background structural details. 

The temporal background noise is generated by the fluctuation in the number of 
incident photons at the detector. This noise usually is referred to as background 
photon noise and is not directly related to the target. Detectors sensitive enough to 
experience background photon fluctuations nolse are called BLIP ( Background 
Limited Photon Detectors ). 

In this work we are interested in the spatial component of the background noise, 
produced by the structural content of the background and mainly related to the target. 
Therefore, from now on whenever we refer to background noise we will mean the 
spatial component. 

Since a target 1s always surrounded by some sort of background it is obvious that 
radiation reaching our detector from this background is confused with the radiation 
emitted or reflected from the target. This 1s a serious problem that complicates the 
detection process and must be well understood for designing a high perforance optical 
or infrared system. [Ref. 1] 

Background radiation is due to thermal enussion and reflected or scattered 
incident radiation. Many parameters affect this kind of radiation. Of high importance 
among them are the changes in emissivity and reflectance of the various background 
surfaces under varied environmental conditions. All the above make the background 
problem complicated and the modeling of a universal background function impossible. 

Background radiance can be considered mainly as two categories: 

1. Terrestrial background 
2. Celestial background 
Terrestrial background refers to radiance seen by infrared systems working within 


the atmosphere of the earth. This kind of background consists of radiation that several 
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surfaces of the earth are emitting or reflecting ( oceans. clouds, vegetation, etc. ). It 
also includes sunlight scattered by particles in the atmosphere and atmospheric 
molecular emission ( sky radiance ). 

Celestial refers to that background seen bv a system working in space and 
consists of radiation from the stars and the sun reflectance from planets and space 
debris. 

Of particular interest in this work are the background radiances of the sea surface 
and the sky. Detailed discussion of these backgrounds and their spectral radiance 
characteristics, especially for the spectral region of 8 to 14 tm, where we have taken 
our measurements, is given in Chapter IIT. 

Since target and background radiances travel through the atmosphere to reach 
the detector, they are heavilv influenced by the atmospheric absorption under various 
meteorological conditions. Detailed analvsis of the atmospheric effects on thermal 


radiation is given in Chapter IV. 


B. BACKGROUND RADIATION CONTRAST 

The target-to-background temperature difference is an important factor i the 
target detection probability of an infrared system. The higher the temperature 
difference between a target and‘its background, the greater the distance and the 
probability of detection. 

This teniperature difference or radiance difference is usuallv expressed with the 


quantity Radiation Contrast, Cr. defined in [Ref. 2] as 


an - Wp 
Cr eae 
Wer + Wp 


where W+, and Wp are the target (T) and background (B) radiant enuttances (W. m-). 
The 8 to I4 ftm wavelength radiation contrast versus target-to-background 
temperature difference for different background temperatures is shown in Fig 2.1. 
When the temperature difference between the target and its background 1s very 
small, or goes to zero, the radiation contrast goes to zero and the detection becomes 
very difficult if not impossible. This effect is known as the washout effect and happens 


for “anibient temperature” targets twice a day due to the solar cycle. 
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Figure 2.1 Radiation contrast for the 8-!4 um region. 


In [Ref. 3] Hudson describes an experiment Witcre tie mmportanceso) tie =eaniaras 
effect is shown. In that experiment a vehicle parked in an open field was viewed by an 
infrared optical system over a period of 24 hours. During that observation. contrast 
Variations were found. In the afternoon the contrast was positive. since the vehicle 
heated by the sun was warmer than tts background. Later in the early evening tours 
the contrast was higher. since the vehicle with larger thermal capacity cooled more 
stowly than the background. As the night progressed. the contrast passed through zero 
and went negative, since the vehicle continued to cool rapidly and at some time became 
colder than its background. During the morning the negative contrast was still present 
because the background was warnung more rapidly. Aiter some hours the feanmesar 
the vehicle was enough to give a period of zero Contrasmamd then timmedsto) positic 
values again. 

Hence, twice in a 24 hour period we are lachineysullicient contrast leveluecruaan 
the target and its background, which means that the detection of anv target at these 
periods is not possible. and since this effect can not be elimmated we have to utilize 


other methods of detection. 


C. BACKGROUND SUPPRESSION 

Target and background have different spectral radiances and different geometries, 
differences that make it easier for them to be distinguished. Whenever a spectral 
radiance difference exists between a target and its background an optical spectral filter 
can be used to reduce the unwanted background signals. Even if spectral. filtering 


Increases the target to background ratio, this is rarely sufficient to render a system 


18 


Operational since usually the irradiance of the background is much higher than that of 
the target. A common characteristic of targets is that thev are alwavs smaller in 
angular extent than their backgrounds, a typical example being a ship against the sea. 
The differences in geometry between targets and background can be expressed as 
differences between their spatial frequency components. For such geometrv differences. 
spatial filtering provides a good suppression of backgrounds of this extent relative to 
the target radiance. 

Spatial filters include all types of reticles discussed in detail in [Ref. 3]. Spatial 
filtering 1s supplemented by spectral filters that define an effective spectral bandpass 
for the system and electrical filters that work on an average background threshold level 
to chop the background de signal. 

[t.is easily understood that spatial filtering is the basic suppression scheme 
mainly giving discrimination of the target in an extended background. But for many 
mulitarv systems. depending on the tasks they perform ( identification. tracking. etc.). 
this is not enough. We have interference from other smaller kinds of background such 
aS Uicited Objects in the Search area, for examiple buildings, rocks, clouds, etc, or 
unexpected objects crossing our detector’s field of view such as birds, different kinds of 
Simms, etc. 

Therefore. more sophisticated and complicated spatial filtering systems have to 
be used for all the different purposes. These svstems are supported by improved 
electronics and computer software, using different implementing techniques to mininuze 


the background effects and optimuze the svstem’s performance. 
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IH. BACKGROUND RADIANCE 


A. MARINE BACKGROUND 
1. Introduction 
The sea surface infrared radiance (L,) is the sum of the radiance caused by its 
thermal enussion (L..(1-R)) from below the surface and the radiance due to reflected 


incident radiation (L;.R) (R: reflectance of the sea surface). 





ie a. Pees) as Deak 


Figure 30) Séasumice pacianice. 


Major factors that determine the Manne aacheromid cliipdcicr sane saa. 
following (Ref. 4]: 
a. Infrared optical properties ol scanvarer 


b. Geometrv of the sea-surface and wave slope distribuuon 


QO 


Sea water temperature distribution 


SS 


Properues of the bottom material . 
2. Infrared Optical Properties of Sea Water 
Thermal radiation from the sea surface occlirs at wavelenuths having tiem 
maxima close to the 10 ttm resiom, AS welcamescomineitoms.2 ton 1 nOrme cl )enen 
these wavelengths Water is essentially (Opaqie im tiienmes smo! eOUiem =): amine 
wavelengths above 3.5 [lm only tie Upperlavers:ot tie seamolsthic hme. ssa ouue Neca 


and less, deternune the sea surface thermal enussion. 
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Figure 3.2 Sea water absorption for selected wavelength regions. 


The radiation emitted from the sea surface is dependent on the underwater 
temperature profile. and is little affected by surface contamination. Since the enussivity 
is equal to the absorptivity the sea surface radiance is dependent on the distribution of 


the absorption coefficient K of the sea water as shown in Fig 3.3. 
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Figure 3.3. Absorption coefficient K of sea water versus wavelength. 
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Subsurface scattering of incident radiation is not significant in the infrared 
range. though it happens in the solar range (0.35 - 3.0 ftm). In that region the 
absorption coefficient has a very low value and radiauon penetrates deeply into the 
water and 1s scattered. 

Reflected incident radiation is composed of solar radiation and sky thermal 
radiation . Solar radiation donunates at short wavelength and of course 1s present only 
in the daytime. In the region from 2 ttm to 4 um reflected solar emussion especially at 
small angles of incidence is greater than the thermal enussion of the sea surface. In 
comparison to this region only about 0.08 percent of the sun radiant energy is found in 
the region of § fim to [4 Jtm wavelengths, which means that only reflected sky thermal 
radiation influences this region by day and night. 

The aniltence ol retlected a auido was recognized by Me Alister [Ref. 6] 
who found. for a clear sky at night, an apparent sea surface temperature a few tenths 
ofa °C colder than the actual temperature. 

Reflectance, transmittance, enussivity and index of refraction for water. are 
shown in Fig 3.4 to Fig 3.7 (Ref. 4]. Flere we have to mention that distilled water and 
sea Water transnuttance do not have any significant difference. 

For the purpose of this thesis. where swe have worked im the $10 1] sam 


wavelength region, We have taken into account only reflected sky thermal radiation. 
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Figure 3.6 Reflection from water surface at different incidence angles. 
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3. Geometry of the Sea Surface 

The slope distribution of the waves is an important factor tn the reflection and 
refraction of incident radiation on a roughened sea surface at various wind speeds. In 
Fig 3.8, taken from [Ref. 7]. is shown the reflectance of solar radiation from a flat sea 
surface (o = 0 ) and from seasurface roughened by a Beaufort 4 wind ( o = 0.2 ), 
where o is the square root of the variance of the wave slope distribution. In this figure 
we can see that for an average rough sea ( wind speed of Beaufort 4 ) the reflectance of 
the sea near the horizon goes to about 20% where as the reflectance of the tlatesumiice 
is about 90°%> and above. Consequently we can see how high the emissivity ( 1 - R ) of 
the roughened sea surface is near the honzon. Only at higher angles of incidence, about 
40°, does the reflectance of the roughened sea surface go slightly higher than the 


reflectance of the flat surface. 
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Figure 3.8 Reflectance froma flat ( o = 0 ) and roughened sea surface (o = 0.2 ). 


In [Ref. 7] Charles Cox and Walter Munk studied this phenomenon for solar - 
radiation in the visible range. They assumed unpolarized sun light reflected once no 
double scattering. They used the index of refraction for uncontaminated sea water as n 
=ieooo. 10rmediculatins the reflection coefficient. Information is lacking for the 
infrared wavelength region and for night conditions. but since the sky radiance 
variation is almost the same during the day and night, and the reflectance of the sea 
water in the visible region of electromagnetic spectrum at zero angle of incidence is 
almost equal to the average reflectance in the infrared region from 2 tm to 15 um ( see 
ce rm mulesceciiycs aresconsidered helpful. For the purpose of this thesis. we have 
used these curves shown in Fig 3.8 from [Ref. 7] after making a correction for the index 
of refraction in the infrared range ( 8 to 14 tm ) we have used mn our work. A detailed 
discussion of this is given in Chapter VIII. 

iiemenO oem itietiiese Curves is that for the dashed portions near the 
horizon, above 80° angle of incidence. the function berng used is not strictly applicable 
because of shadowing and multiple reflections. Also the lower and upper branches of 
Pew woe — Um seclive tepresent two assuniptionis used im [Ref. 7] for the multiple 
fceleeionmeiieci. [ne basic fimetion used works for light reflected only once as we 
mentioned before. True values are expected according to C.Cox and W.Munk to lie 
between the indicated limits. 

Miieweiives im Fie 3-9 show the sea suniace radiance (jt) divided by the sky 
radiance \.(0) at the zenith, due to sun light scattering by air molecules, as a function 
of the vertical angle jt. for a clear skv and an overcast skv, and 1s also taken from 
ict. 7). 

4. Sea Water Temperature Distribution 

The sea surface temperature determines the enutted thermal radiation of the 
sea. This temperature is about 0° C (mainly about 4°C) in arctic regions and increases 
reaching about 29°C near the equator. Currents may produce anomalies of some 
degree in temperature. However the thermal emission of the sea surface has its 
maxinium peak near the 10 pm wavelength. As we have said before mn the infrared 
optical properties section, the temperature of the upper lavers of the sea surface. of 
thickness about 0.01 cm. characterizes the thermal enuttance. 

Under evaporation conditions the teniperature of the upper layer of 0.1 mm of 
the sea surface was measured [Ref. 4] as colder than the water temperature a few 


~ 
e 


centimeters below the surface of about 0.6°C. In Fig 3.10 [Ref. 4] are shown 
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Figure 3.9 “Sea suriace radiance diniicd by shar oiicesacezemiune 


measurements of typical conditions showing the sharpest gradient in the upper layer of 


aan 
>. Properties of the Bottom Materials 
These properties are important in regions of short wavelengths where 


radiation can penetrate deep into the water causing subsurface scattering. Such 


properties are nol important for the work done im this (hesismsinee: we wiaemn och comm 


the region of longer wavelength ( 8 tm to [4 bm ). 
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Figure 3.10 Thermal structure of sea boundary layer. 


Peace eeyRAL RADIANCE OF THE SKY 
1. Introduction 

Skv background radiance in the infrared is caused by uvo different factors, the 
scattering of radiation from the sun and the thermal enussion of the atmospheric 
constituents. The scattering factor 1s important in the short wavelength range. close to 
fier isiOlewmand IS presen: only im daytime, The thermal emission happens at 
Pewee cne.ns 1oneer than 4 ini and is present dav and nicht. 

An idealized form of this contribution to the sky spectral radiance ts given in 
Fig 3.11 (Ref. 8]. This figure shows the spectral radiance of a blackbody at 300° K, 
close to the temperature of the earth, and the spectral radiance of a blackbody at 
6000° K. close to the temperature of the sun. On estimating the sky radiance one can 
assuine that the radiation from the sun is diffusely scattered and henee the sky appears 
to have a uniform radiant energy. This kind of maximum brightness, which has an 
average radiance about 2 107° times the radiance of the sun. can De ASELONIMiaIed OF 


a bright cloud shown as the lower dashed curve in the figure. But normally the sky 


does not scatter this large amount of radiation; a nominal value is about one-tenth of 
this maximum spectral radiance, shown as the solid curve on the left, named “clear 
sky". So we can describe the combined effect of the thermal radiation emitted at the 


temperature of the earth and the scattered sunlight as combination of the two solid 


curves shown in this figure. 
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Figure 3.11 An idealized spectral radiance of the sky. 


2. Atmospheric Emission 
a. General 

The atmospheric enussivity in the infrared depends on the amount of water 
vapor, carbon dioxide and other gases which give rise to absorption bands due to their 
vibrational and rotational transitions. Since the absorptivity equals the emissivity 
( Kirchhoff's law ), we can predict the atmospheric emission from the absorption 
Spectra of these gases. 

b. Spectral Radiance of Clear Sky 

Skv temperature and elevation angle are the main factors determining the 
radiance of a clear sky. The elevation angle determines the length of the atmospheric 
path and hence the emissivity, while the temperature determines the blackbodv 
enussion characteristics. 

The effect of the elevatron angle on spectral radiance is presented in Fig 
3.12 [Ref. 8]. At low angles of elevation ( near the horizon ) we can see nearlv 
blackbody radiation at the ambient temperature. But as the elevation tncreases above 
the horizon and the atmospheric path length decreases the optical density decreases 
and the enussivity 1s no longer blackbody ( € = 1 ) everywhere. In the regions 
centered about 6.3 [fm ( water vapor absorption band ) and about 15 um ( carbon 
dioxide absorption band ) enussivitv is very high, and the skv remains blackbody, even 
at the zenith ( 90°). But between these regions, in the region from about 8 jim to 13 
jin the enussivitv is poor and the sky is no longer black. The general background 
emussion in this region 1s due to the wings of the water vapor and carbon dioxide 
bands. The peak that we see at about 9.6 fim is due to ozone enussion, which 
Increases as the elevation angle decreases from the zenith. 

Fig 3.12, taken from [Ref. 8], shows spectra of clear sky radiance measured 
in September, at night. from a height of 11000 feet above sea level. at the ambient 
temperature of 8°C. The angles of elevation in sequence are 0°, 1.8°, 3.6°. 7.2”. 
14.2°, 30° and 90°. 

The effect of the ambient temperature is very strong for this kind of 
radiation, as shown in Fig 3.13 . This figure is also taken from [Ref. 8] and shows the 
zenith sky spectral radiance for two different ambient temperatures. The dashed curves 
are the spectral blackbody radiances corresponding to the ambient temperatures for the 
two separate measurements at 27.5°C and 2.5°C. These measurements were taken with 


different amounts of air mass, which affected the emussivitvy and hence the radiance of 
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Figure 3.12 Sky radiance for several elevation angles above horizon. 


the spectrum in the regions with poor emittance [Ref. 3]. The low temperatane 
measurement was taken at 14110 feet above sea level and the higher temperature 
measurement at 6000 feet above sea level. The negative values shown are due to 
measurement error. 
c. Spectral Radiance of an Overcast Sky 

Overcast skies have been found to be very good blackbody radiators ( € = 
|), having a temperature at low altitudes within one or two °C from the ambient 
surface temperature Reiesi 

As mentioned before, water in the infrared is a good absorber (see Fig 3.3) 
with a small reflectivitv. Clouds mostlv at lower altitudes consist of water droplets with 


radii from 1 jtm to about 10 [tm at concentrations above 300 cm” with temperatures 
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Figure 3.13 Variation of zenith sky radiance with ambient temperature. 


not far from the surface tenyperature [Ref. 9]. So infrared radiation incident on a cloud 
water droplet 1s highly absorbed and partially reflected and transmitted. Hlowever due 


to multiple scattering from all neighbouring droplets in the cloud the transmittance of 


Becee @ucNeecs eSsentiall, to zero. 
Flenee we can sav that clouds in the injrared region are good absorbers and 


consequently good emutters, causing overcust skies to behave as good blackbody 


FACtalOnrs. 


d. Spectral Radiance of Clouds 
Clouds with thickness greater than 100 m and clouds with water droplet 
; : 4 ) 
( or ice crystal for higher altitudes ) counts greater than 0.5 g/cm~ do not transmit. 


Thus absorptivity, and consequently enussivity, is high, approaching one, which means 


that clouds show a verv good blackbody behavior {Ref. 10}. 

Fig 3.14 [Ref 10] shows the enussivity of cloud lavers versus wavelength 
and thickness AZ. Also in Fig 3.15 [Ref 10] is shown how the integral enussivity (CL), 
transnuttance (1) and absorptance (A) changes with thickness AZ. Here we can clearly 
notice that at thickness AZ above 60 m the transnuttance goes to zero and above 100 
mi the enussivitv goes to zero. 

However, looking at the spectrum of the under side of a cumulus cloud in a 
clear sky ( Fig 3.16 ) at an elevation angle 14.5° [Ref. 8] we do not observe a 
blackbody spectrum . Since the emissivitv of the cloud 1s expected to be high and since 
the sky is partiallv transparent in the region of 8 - 13,tm for the elevation angle of 
14.5°, as shown before, we might expect to see clearly the emission of the cloud in that 
spectral region. From the figure we see that the emussion in this region is characteristic 
of a blackbody at -10°C which is reasonable for the temperature of the cloud ( lower 
dashed line in figure ). What is really happening in this figure is that we have a 
combination of the atmospheric enussion at the 6 and 15 Jlm wavelength regions from 
the path below the cloud together with the -10° blackbody enuission from the cloud in 
the wavelength region at $- 13 tm. 

The effects of thin cirrus clouds in a clear sky modified situation is 
presented in Fig 3.18. Here we can see that the emission in region 8 -13 jtm has been 
raised in coniparison with the clear skv situation ( see Fig 3.13 ) , and also the ozone 
enussion 1s reduced since it comes from a region higher than the clouds and is 
attenuated passing through them. 

Fig 3.17, from [Ref. 8], shows spectra mieastired in June at sea lewe] wank 
an ambient temperature about 28°C. The elevation angles are 0°, 18°, 3.6 tae 
14.5°, 30° and 90°. 

3. Scattering of the Sky Radiation 
The position of the sun relative to the area of the observed sky and the length 
of the scattering path due to the elevation angle affect strongly the amount of scattered 
radiation. We have a large amount of scattering close to the sun direction as Well as an 


increase in scattered radiance with a decreasing elevation angle. Decreasing the 
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Figure 3.14 Spectral emissivity of cloud lavers of thickness AZ. 
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Figure 3.15 Integral enussivity (E).absorptvity (A).transmuttance (T). 


Pewamomeanete we have a longer atmospheric path from: which scattering can be 
produced. But verv close to the horizon we have a reduction in the scattered radiance 
mainly because of water vapor absorption which reduces the effective observation path. 


This kind of radiation 1s not important for the work done in this thesis since 1t 


occurs at short wavelengths. 
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Figure 5.16 Radiance of the under side of a cumulus cloud. 
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Figure 3.17 Spectral radiance of sky covered with cirrus clouds. 


IV. ATMOSPHERIC TRANSMITTANCE 


A. ATMOSPHERIC PROPAGATION 

Thermal radiant flux from a target has to pass through the atmosphere of the 
earth before it reaches an infrared sensor. During this passage, the thermal radiation of 
the target 1s attenuated because of selective absorption due to several atmospheric 
gases and scattering by small particles suspended in the atmosphere.such as aerosols. 
rain, snow, smoke, fog and haze. Both absorption and scattering depend on wavelength 
(ei); 

In general we call atmospheric extinction the process by which the radiant flux 
from a target 1s attenuated passing through the atmosphere. This can be expressed as: 

T (A) = e H(A)R | 

where TA(A) is the atmospheric transmission, 6(/) is the extinction coefficient and R is 
the range. The extinction coefficient is the sum of the absorption coefficient a(?.) and 


the scattering coefficient 7(A) 


O(A) = ah.) + YA). 

Both absorption and scattering cefficients are sums of molecular and aerosol 
components [Ref. 2]. [In the infrared region the absorption process is of greater 
importance (particularly in the 8 - 14 jtm ) than the scattering process which 1s 
important in the visible and near infrared region. 

Water vapor, carbon dioxide and ozone are the most significant absorbing gases. 
having broad absorbing bands centered at the specified wavelengths: water vapor ( 2.7. 
eeoe wii wearvon dioxide ( 2.7, 4.5, 13 iim ). ozone (4.8. 5.2, 6.3 ftm ). The bands 
of carbon dioxide at 2.7 and 15 jtm and water vapor at 6.5 ttm limit the atmospheric 
transmission in the infrared to two atmospheric windows at 3.5 to 5 tm and at 8 to l4 
jim. 

Infrared systems are designed to operate through these windows, and it 1s very 
important for the performance of our system to have a prediction of the atmospheric 
attenuation for the wavelength and the environmental conditions in which we are 
Operating. This prediction process is a complicated problem and has been faced by 


computer modeling codes. 


toes 
Can 


Such models are the high resolution models ( HI TRAN) and the low resolution 
model ( LOWTRAN). The high resolution model codes provide direct computation 
with high resolution single frequency results. They are used to predict atmospheric 
attenuation for laser radiation. The low resolution model codes are used for a broad 
band prediction where the thermial svstems work. 

LOWTRAN code was used to predict atmospheric transmittance for the work 


done in this thesis and 1s considered in the next sections of this chapter. 


B.  LOWTRAN GgCOpE 

1. General Description 
LOWTRAN 6, the latest edition of the LOWTRAN simulation codes: is a 
FORTRAN computer curve-fitting program, relating empirical precomputed 
attenuation data to simple meteorological and environmental inserted parameters. It 
was designed to calculate average atmospheric transmittance and thermal radiance 
over a specified broad band, ranging from 359 to 40000 cm’! ( 0.23 to 28.57 {ni ) wath 


averaging bandwidth in steps of 3 em‘!, 


a resolution of 20 cm” 

The average transnuttance 1s a function of the total extinction, the product of 
the transmuttances due to water vapor line absorption, water vapor continuum 
absorption, uniformly nuxed gases line absorption, nitrogen continuum absorption, 
aerosol absorption. aerosol scattering and molecular scattering. A consequence of the 
averaging over a finite spectral band is that the rapid fluctuations of the atmospheric 
transmittance in that band vanish. 

Generally we can sav that the LOWTRAN code is a fast but not very 
accurate computer simulation code. 

2. LOWTRAN 6 Input Data 

The LOWTRAN © code 1s activated by submitting four cards of sequential 
data. Each card provides data to be selected and asks for data values to be inserted as 
follows: 

The first card provides a selection of six atmospheric models. These models 
are the tropical atmosphere, midlatitude sumuner, nudlatitude winter, subarctic summer, 
subarctic Winter and the 1962 U.S. standard atmosphere. Also in this card we have to 
indicate the tvpe of atmospheric path ( horizontal, vertical or slant ), and determine the 


execution mode ( transnuttance - radiance ). 


(3 
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The second card provides a selection of the type of extinction for nine 
bountary- layer aerosol models and gives a default meteorological range. These models 


are the rural (vis = 23 km). rural (vis = 5 km), navy maritime (computes visibility). 


Cay 


maritime (vis = 23 km). urban (vis = 5 km), tropospheric (vis = 30 km). advection 
fog (vis = 0.2 km), radiation fog (vis = 0.5 km) and a user defined model (vis = 23 
km). Also here we have to specify the air mass character for the navy maritime model 
and some meteorological conditions. 

The third card defines the geometrical path and we have to insert some 
parameters lke initial and final altitude, initial zenith angle and length of our path. 

In the fourth card we have to specify the spectral range over which we want 
data for atmospheric transmittance or radiance and the frequency increment over 
which we want our data to be averaged. 

A detailed description of the functions of these cards and the required input 


data is given in [Ref. 11]. 


C. LOWTRAN 6 CODE USED IN THIS WORK 

For the work done im this thesis the LOWTRAN 6 code was used to find the 
atmospheric transmittance for the different meteorological conditions we experienced . 

For our experimental conditions, mostly overcast skies and ambient temperatures 
about 10°C to 16°C we selected the midlatitude winter ( 45° N, January ) atmospheric 
model for horizontal path. Also for the purpose of this work the navy maritime 
bountary layer aerosol model was selected. And thus because this model takes tmto 
account the existence of three components in deternuning the aerosol population. First 
is the continental component which is sensitive to the amount of time the air mass has 
Spent over the sea; Second is the stationary component which depends on the wind 
history and represents aerosols formed bv high winds and whitecap conditions, and the 
third one is the fresh component which is characterised by the current wind speed and 
consists of whitecap produced droplets. For our geometrical path we mserted target 
altitude 10 m and horizontal path length 650 m. Our spectral range was 710 - 1250 
cm’! (8 - 14 pm ) with frequency increment 5 cm’). 


The results of these calculations are discussed in Chapter VII. 


oF 


V. FAMILIARIZATION WITH THE EQUIPMENT 


Aen GENERAL 

As mentioned in the introduction. all our data for this work were collected, 
stored and processed using the AGA Thermovision 780 svstem, supplemented with a 
digital microcomputer. 

The system is designed for both thermal imaging and image analysis. It 
combines real infrared scanning with thermal image measurement capabilities and 1s 
supplemented by a microcomputer for digital image recording, storing and processing. 

The particular system we used was a dual scanning svstem operating in the 3 - 
5.6 fim and.8 - 14 jtm windows with two different detectors. For the purpose of this 


experimental work only the $ - 14 [im scanner was used. 


B. AGA THERMOVISION 780 
1. Basic Description 

The AGA Thermovision 780 system consists of a real time scanner, a black 
and white monitor chassis display and a digital microcomputer. 

The real time infrared scanner is the unit that converts the thermal radiant 
energy emitted from a target to electronic video signals. These signals are amplified 
and transmitted through a cable to the display monitor chassis where with further 
amplification thev produce a thermal image on a black and white screen. The amplified 
signal can also be fed to a color monitor through an interconnecting device in order to 
have a real time thermal image presentation. Finally video and data signals are 
transmitted through an interconnecting link to a microcomputer. Here the thermal 
image is desplaved in color and can be recorded, stored on a data disk and processed 
for temperature evaluation using color codes. 

2. The Infrared Scanner 

There are two Kinds of scanning unit used in this svstem, the short wave (SW) 
scanner that operates in the 3 to 5.6 Jtm spectral band and the long wave (LW) scanner 
that operates in the 8 to 14 jtm band. Scanner systems are designed with single or dual 
scanner units which cover both bands simultaneously. Our scanner as we have already 
mentioned was a dual tvpe. Such a scajimer svstein is Showa in Pic Sal: 

Each scanner unit consists of the following parts : 


a. Electro-optical scanning mechanism 
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Infrared detector 
Liquid nitrogen dewar for cooling the detector 
Control electronics and preamplifier 


A simplified block diagram of the scanner unit is shown in Fig 5.2. 





Figure 5.1 Dual scanner svstem. 
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Figure 5.2 Simplified scanner unit block diagram. 


Electromagnetic thermal energy emitted from a target 1s focussed by infrared 
lenses into a vertical prism. The lenses for the short wavelength band (SW) are made of 
silicon and those of the long wavelength band (LW) are made of germanium. Both 
kinds are covered with proper antireflecting coatings. The vertical prism is rotated at 
[SO rpm and its optical output is passed through a horizontal prism which rotates at 
[S000 rpm. 

Rotating slotted discs control the rotation of both prisms and are 
electronically connected to the horizontal and vertical triggering circuits to provide 
horizontal and vertical triggering pulses for the monitor. The horizontal and vertical 
prisms are synchronized in ‘such a way that four fields each having 100 horizontal 
scanning lines produce one interlaced frame. Of these lines only 70 hnes per field or 
280 per frame are used as active imaging lines. 

The output from the horizontal prismi passes through a set of relay optics 
containing a selectable aperture turret with eight aperture discs between [1.8 and f 20 
(f/1.8 dise was selected fer Glip Case ait (hen OassceathrOuUgmeaiitemmmit ( no filter [on 
our case ) and finally is focussed on a single detector located on the wall of the dewar 
chamiber. The detectors are a photovoltaic Indium-Antimonide (InSb) for the short 
wavelength band and a photoconductive Mercury - Cadmium - Telluride (HgCdTe) for 
the long wavelength band. Both detectors are cooled with liquid Nitrogen which 
niaintains the chamber at a temperature of -196°C. | 

The signal output from the detector is preaniphfied within the scanner unit 
and the video signal produced is applied through a-cable to the black and white 
monitor chassis [Ref. 12]. 

The arrangement of the electro-optical comiponents of a scanner unit is shown 
in Pic Seater 

3. Black and White Monitor Chassis 

The video output signal from the scanner unit is apphed to the input 
connector of the black and white chassis. This signal after a further amphfication and 
processing within the monitor is apphed to the display screen. 

Outputs from this monitor chassis are used for a color monitor connection 
and the interconnecting link for the nucrocomputer connection with the system. 


A Black and White monitor chassis is shown in Fig 5.4 [Ref. 12]. 
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Figure 5.3 Arrangement of electro-optical components in a scanner. 





Figure 5.4 Black and White monitor chassis. 


4. Technical Data 
Technical data for the AGA Thermovision 780 were found in [Ref. 13] and are 
presented in fapieue 

5. Digital Image Processing System 

This system supplements the AGA Thermovision 780 svstem and consists of 
an interconnecting link and the microcomputer which utilizes the systems software 
DISCOe 

a. Interconnecting Link 

The interconnecting link is a digital on-hne interface for direct interfacing 
between the Thermovision 780 and the IF 800 model 20 microcomputer. The main 
purpose is to digitize the image, hold the image data during the time required for 
transfer from the Thermovision to the microcomputer, and control the data 
_transferring speed. 

b. LF 800 Nlicrocomputer 

IF SOO is a business level personal computer which uses a Z-SOA . 4 MHz 
clock mucroprocessor for its CPU. It comprises a keyboard, a printer. a display unit 
and a muni floppy disk drive. 

6. DISCO 3.0 Software Description 

The software used bv this system consists of three main programs the DISCO, 
the IMAGE PROCESSING and the UTILITY programs. Various subprograms as 
shown in the DISCO 3.0 block diagram Fig 5.5 [Ref. 14] supplement the main 
programs. 

The main program DISCO was the one mostly used for this work to record, 
store and process our data. Using this program infrared images were transferred from 
the Thermovision 780 svstem to the nucrocomputer. The subprogram Fl Record 
[mage was used to record these thermal images and displav them in color on the 
computer screen. Also this subprograni was used to enter or alter our parameters 
(emissivity, atmospheric transmittamee etc.), to Select and enter the Scanner 
identification characteristics used and store the images on a data disk. The subprogram 
F2 Evaluate Irage was used to evaluate the temperature measurement of the image 
using color scaling with tables of 3 to 36 different color combinations. Each color was 
assigned a particular temperature range. The standard table of the svstem uses $ colors 
and was the one used in our case. This subprogran. was also used to displav 


histograms with a frequency distribution of the temperature values, to display single 


TABLE 1 


aoe HERRON ISTON 780 


FREFORMATCE 


lee NICAL CHARACTERISTICS 


Spectral Range : 3 - $.6 tm and/or 8 - 14 jm 
Frame Rate : 6.25 sec”! 
Field Rate : 25 sec} 
Interlace : 4: 1 
REFLACEABLE FORE OFTICS 
FOV Azimuth * t Je 
FOV Elevation *  s: Je 
IFOV) Azimuth : 1.1 mrad 
IFOV Elevation : 1.1 mrad 
NEAT OL lZ%e Mate zz2cc 
MPT t < "Or rec 
Dynamic Range : -20 to 900e¢c¢ 
OFTICAL DATA 
Effective Aperture Area 24 cm? 
Aperture Diameter So em 
“Effective Focal Length : 9.9 cm 
f/number : ae 
DETECTOR 
Type SW Fhotovoltaic InSb 
LW : Fhotoconductive HgCdTe 


Humber of Elements : 1 


Feak Wavelength SW: 5 


LW : i 


COOLING SYSTEM 


Lquid Nitrogen 


* 


Hote: 


FOV Seagn0. 5, 1X7, 


jim 


QO jtm 


-196e¢c¢ 


The system can take five different fore lenses with 


12X12, 20%20, and 40X40. The lens 


we used was the one with FOV 7X7. 
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isotherms, magnify the image, transfer a new image on the display from those which 
had been stored on the data disk and print data such as temperature profiles using the 
system printer. For the last case our system had some limitations since the svstem’s 


maximum capabilities are available only on a color printer which was not available. 


C. THERMAL MEASUREMENT TECHNIQUES 
1. Introduction 

As we know the radiant flux density received by an infrared svstem does not 
have a linear relationship to the object temperature. It is always mixed with reflected or 
enutted background radiation and is attenuated Dy Ine atmi@sedere. [herelore. certain 
calibration and correction procedures have to be applied to a particular measurement 
in order to have an accurate teniperature evaluation. 

The AGA Thermovision 780 system measures thermal radiation between 
Spectral limits as mentioned before. The numerical value of the received radiation is 
called the THERMAL VALUE and-‘is measured in /sotherm Unirs (1U) which is an 
arbitrarv unit svstem. This thermal value is proportional to the received photon flux 
enc DUlsicmet Piopermona: tol the temperature of the target. Therefore. we need a 
calibration function as a graphical curve or a computer program to relate the thermal 
value ( isotherm units ) with the teniperature of the target. 

2. Calibration Curves 

Cahbration curves for all aperture sizes used bv the system in the long 
wavelength ( LW ) spectral band, in which we worked, are shown in Figs 5.6 through 
5.8 [Ref. 12]. These curves have been obtained assuniung enussivity € = 1 (blackbody 


radiator ) and neglecting atmospheric damping T 1. For situations different from 


w 
ie abOve if isObs1a@us that we have to apply correction techniques. 
These calibration curves are accurately described bv the following equation 
fake [rom | ers 2] 
A 
{= —————__ (egnony 
celBI)_| 


~ t 


where I is the thermal value in isotherm units ( IU ) for the absolute temperature (T) 
in °K. A,B,C are calibration constants that depend on the aperture. filter. scanner etc. 

This equation was used by the microcomputer program to calculate the 
temperature of the target for a known thermal value after the calibration constants 
PeEe tor tnese siown in Table 2 had been inserted. 
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Figure 5.6 AGA LW calibration curves from -20°C to 30°C. 
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Figure 5.7 AGA LW calibration curves from -20°C to 150°C. 
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3. Basic Thermal NMleasurements 

There are two basic methods we can use to measure temperature using the 
thermoviston 780. The direcr, and the relative nveasurement method. 

In the direct method of measurement the system built in capabilities (clamping 
and temperature compensation) permit measurement of the temperature of an object 
Without using any other reference source of radiation. 

In the relative measurement method the radiation of an external reference 
source with known temperature and enussivitv 1s compared with the radiation received 
from the object to obtain the object temperature. 

The direct method 1s easier to use but not as accurate as the@relative method; 
which, however, has the disadvantage that 11 neeqs=an accurate Trelehemec sOulrcetoune 
available. In our case the direct method was used since we didn’t have anv reference 
source With known teniperature and enussivity located near the target. 

4. Direct Measurernent Procedure 

In this subsection we will describe briefly the proceedure for the direet 
measurement method, using the Thermovision first alone and then supplemented with 
a microcomputer as was done to collect data for this work. 

To aid tn understanding this measurement procedure, the monitor chassis 
front panel controls and indicators are presented in Fig 5.9 [Ref. 12]. Subsequently a 
short description of the operation and function of the controls and indicators referred 
to in this measurement procedure ts given. 

a. THERMAL RANGE (item 9 ):. This is a nine position switch calibrated 
between 2 and 1000 isotherm units. This switch selects the thermal span of 
interest. 

b. THERMAL LEVEL (item 10) Whis control sets tine thermal eve: of thiemimac] 

c PICTURE MODE “NORMIATS 3 itenii2om melieselects tie 110rn ice 
Cea Where the cooler object temperature appears as black and the hotter as 

d. ISOTHERM LEVEL 1 (iteny 135)) Wis contrelstuins on themsetherm funetion 
and identifies the discrete thermal levels on the target. The isotherm funetion 1s 
indicated by saturated white areas. which Tepresent areas of the Same 
teniperature. A marker on the ventical isetienmi scale indicates the relate 
position of the selected level. 

e. ISOTHERM SCALE (item 16 ): This is a scale calibrated from ve to +0558 


reading of the isotherm marker multiplied by the thermal range gives the 
relative thermal value tn isotherm units. 


aw, 


For both cases either with the thermovision itself or supplemented by the 
mucrocomputer, starting with picture mode in NORMAL, we have to adjust the 
THERMAL RANGE and LEVEL controls to obtain a good thermal picture. The 


fan 7. 9 i 12.10 





dg POWER 8. CONTRAST 

Ze Les 3a. CONTRAST PRESET 
3: FOCUS = THERMAL RANGE 

a) LINE/ BATT 10. THERMAL LEVEL 

2. POWER LED ti THERMAL LEVEL AOu. 
6. PHOTO 2 PICTURE MODE 

2. BRIGHTNESS uc ISOTHERM LEVEL 1 
7a. SRIGHTNESS PRESET 14. ISOTHERM LEVEL 2 


er ISOTHERM WIOTH 
16. ISOTHERM SCALE 


Figure 5.9 Mfonitor chassis front panel controls and indicators. 


picture was shown on the Black/White monitor and on the color monitor used. The 
measurement procedure is shown schematically in Fig 5.10 (Ref. 12]. First with the 
THERMAL RANGE control switch we select the thermal span of interest in isotherm 
umits ( 2 to 1000 ), about the thermal level of the image. In our case thermal ranges of 
a2 ame, 10 isotherm units were selected. Then with the THERMPAL LEVEL control 
Knob we set the thermal level of the image. At the point at which a satisfactorv 
thermal picture is obtained the corresponding thermal value “L” of the image in 
isotherm units, is measured (i1.e., the thermal value of the isotherm scale at zero point. 
Fig 5.10). 

omosrnie thermovision alone. the control ISOTHERM LEVEL | has to be 
adjusted to brighten up the point of interest on the target in view. A marker on the 
vertical scale indicates the relative position of the selected level. The reading of the 
isotherm marker multiplied by the thermal range gives the relative thermal value “1” in 


isotherm units. shown in Fig 5.10. Adding the two values “L” and “i” we have the 
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Figure 5.10 Direct measurement method. 


measured thermal value I’, for the particular point of interest. Under the assumed 
simplified conditions, having tT, = land eé = I, the calibration curve can be used, as 
in Fig 5.10, to translate the measured thermal value I’, to temperature. 

Using the microcomputer as was done to collect data for this work there was 
no need to go through the isotherm procedure. The computer does not use isotherms 
but discrete thermal values. Therefore, we had the isotherm marker at zero. The 
subprogram Fl Record Image was used to record not only the thermal image picture 
but also the corresponding thermal value and thermal range. These values were used by 
the subprogram F2 Evaluate Image to evaluate and present on the display the 
temperature distribution of the recorded image in color scaling, using the standard 


eight color combination table of the system. 
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D. COMPLEX THERMAL MEASUREMENTS 
1. Introduction 

As mentioned before, the AGA Thermovision 780 system considers an ideal 
situation for the values it measures. But as we know there are several factors that 
influence these measurements. So the true object measurement must be derived bv a | 
calculation in which we must take into account alt the factors that influence the 
received thermal radiation. 

2. Variables Affecting the Measurement 
The most important factors that affect the measured thermal radiation are the 
following: 
a. Emussivitv of the target 
b. Background surrounding the target 
Target opacity 
d. Target size 
e. Atmospheric effects 

Most of the objects in our environment are not blackbodies and hence their 
enussivities are less than one. For our case the emissivitv of the R)V “POINT SUR’ 
was found to be about 0.95 as will be explained later. The influence of the background 
was discussed in detail in Chapters I] and III and the atmospheric effects were 
discussed in Chapter IV. 

The target opacity has to be taken into account because some targets are 
more or less transparent to infrared radiation ( glass, plastics, clouds, ete. ). Radiation 
received from objects like these includes radiation transmitted from the background. 
PcOme@emeniissivity oi the target is affected and is hard to estimate. For our case the 
target was considered as an opaque Lambertian surface. — 

Wiiesi@eso! the target is another factorao,be taken into account. Since the 
detector subtends a specific solid angle, if the target is not big enough to cover this 
angle the detector receives radiation from the target background. This will tend to 
affect the temperature difference indicated between the object and the background. In 
our case the target was not big enough to cover the field of view of a 79X70 Iens ata 
distance 650 m. For most measurements it occupied 1.3 to 2/5 of the display screen. 

3. Exact Measurement Formula Derivation 

the thermal radiation that reaches the detector of our svstem is the sum of 

the target-enutted thermal radiation, the radiation reflected by the target and the 


radiation emutted from the atmosphere as shown tn Fig 5.11. 


oe) 


OBJECT 





ATMOSPHERE 


SURROUNDINGS 


Figure 5.11 Radiation conditions in complex mieasurement. 


If we assume an opaque object at temperature (T.,) and enussivity ( €, ), this 


object emits photon flux €)5,. S, is the photon flux of a blackbody at the same ( T,) 


O 
temperature. Since the above flux 1s attenuated by atmospheric transnuttance ( t, ) the 


object radiation reaching our detector is €,Ty4S,. 


Also assuming that all surrounding surfaces are of the same ambient 
temperature ( T,, ) which determunes that the enmussivity is close to one ( €&, = 1 ) we 


can say that the radiation reflected by the object is p,Sy. S, is the photon flux of the 


ad 
surroundings emutting as a blackbody ( €, = | ) at temperature ( T,, ) and p, is the 
reflectance of the object. From the first assumption that the Object 1s opaque we ger 
Pp, = 1 - €, and since the above radiation is attenuated by the atmosphere the 
reflected radiation reaching our detector ts ty( 1 - £4 )Sg. 


Kg Sar where aa 1S 


the photon flux due to atmospheric enussion and scattering. Hence we can write that 


Finally the atmospheric radiation ts given by ( 1 -T 


the photon flux S°, reaching the detector ts 


So = =6'u > en ee ea (ecmes. 2). 


OQ atm 


rd 


This photon flux relation can be converted into a thermal units relation since, 
as we mentioned before, the thermal value (JL ) meastred@ar= tie, 47 emo ision 


780 is proportional to the photon flux reaching the detector itiere! orem vemcame wre 


the relation [ = C X-S where I is the thermal value. C the proportionality constant 
emcees the plioton flux. 

Substituting S = I,C in Equation 5.1 we have the same equation in thermal 
values 


oimetice ee fai! - ey ae i Tee Pann 


(ean O25): 

This equation can be solved for I, to find the exact temperature through the 
calibration curves or can be solved to find the emissivitv E, if we know the exact 
temperature of the object and consequently I, from the calibration curves. Therefore 


POmeediaton 5.5 we have 


V 


O I l l 
I. =-——-(-—-1l1)l,-—(—-I1 JI. , (eqn 5.4) or 
O a atm 
Ao Mean 5 ee 
Yo > tglg - Cl - tylatm 
2 (eqn 5.6) 
emis) 


For our case the computer program calculated for us the exact temperature 
using Equation 5.3 after we had inserted certain values for the parameters enussivitv. 


atmospheric transmittance, ambient temperature and temperature of the atmosphere. 


VI. EMPIRICAL CALIBRATION OF AGA THERMOVISION 780 


A. CALIBRATION PROBLEM 

Trving to process the thermal images we had already recorded. we realized that 
the temperatures the computer was calculating were much higher than expected for the 
particular scene. 

It was obvious that the Thermovision 780 was measuring thermal values ( IU ) 
higher than the values the calibration curves were giving for certain target 
temperatures. 

To solve the problem we had to derive mew calibratiompeuniesstomse sionuiue 
eXact temperature evaluation of the stored thermal images. “At taarerimee nad 
difficulty deriving new calibration constants to be inserted in the computer prograin. 
Therefore, we tried to find a relation between the exact temiperdiuie ro! aeoimecsom 
and the temperature the computer program was calculating using the faulty calibration 
curves for the blackbody thermal values measurements. In this way we had an 
empirical relation as a curve function, accurate enough to correct the temperatures the 
computer program was calculating using the faulty calibration curve on processing the 


thermal images we had stored. 


B. THE CALIBRATION PROCESS 
1. General 
For the calibration procedure we required a blackbody working at about thie 
temperatures we had when the images were recorded. lie air teinperattinesat tlie 
period was 9.7 - l).3°C, the sea striice tempenitiine v75 slike? = ele 3 © Sane 
laboratory ambient temperature, where the calibration procedure was carried out, was 
always 18 - 19° C. Therefore we needed to have a blackbody source to work below the 
ambient laboratory temperature in order to get the calibration relation we wanted. 
2. Blackbody Source 
For the purpose of this work we built our own blackbody source to use in our 
calibration procedure. As we have already mentioned we wanted this source to work 
at about the expected target temperatures, which were below the laboratory ambient 


[emMperaicune: 


SO 


The core of this blackbody source was a cylinder 28 cm long and 7.5 cm wide 
made of thick (0.5 cm ) aluminum. One end of this cylinder was cut at an angle about 
60¢ from horizontal and it was covered with a plate of the same thickness of aluminum 
for nunimizing reflected radiation escaping out of the cavity. The other end was cut at 
90° and it was covered with a same thickness aluminum plate which had a hole of 2 
cni diameter. The interior was painted with mat high emissivity black paint. Using a 
nichrome wire wound unifornuy around the cylinder we could heat our source and by 
pouring liquid nitrogen inside the cylinder we could cool it below ambient temperature. 
The whole source was covered with a kind of fiber glass insulating material and four 
thermocouples were used to check the temperature stability of the core cvlinder and the 
edge plates. [n Fig 6.1 a schematic representation of this blackbody source is shown. 

This source was well tested showing very good blackbody behavior, stabihzing 


and maintaining the temperature We needed quite fast. 





-—7.5cecm= ——— — 23m = 


Figure 6.1 Blackbody source used for calibration procedure. 


5. Calibration 
We wished to find a relation between the exact temperature of an object and 
the temperature calculated by the computer program using the wrong calibration 


curves for the isothermal values measured by the Thermovision 780. 


For this reason our blackbody was set to certain temperatures about the 
expected target temperature at the time of the experiment. After these temperatures 
Were stabilized a measurement with the Thermovision 780 was taken from a distance of 
about 2 m. The thermal value was transferred to the computer and the thermal image 
was recorded, to be processed for temperature evaluation. For the distance of 2 m the 
atmospheric attenuation was neglected ( t, = 1). The above procedure was repeated 
for all three thermal ranges { 2. 5, 10 ) used to collect data. 

The above calibration temperature measurements for the three different 
thermal ranges are shown in the following Tables 3.4.5. Looking at these 
measurements we can notice that there is a linear relation between the blackbody exact 
temperature and the temperature the computer was calculating. Using the method of a 
least squares linear regression we found the equation of the best fitting straight line. 
This line was the calibration curve function we wanted to correct the temperatures of 
the stored thermal images. 

Plotted calibration data measurements and the best fitttng straight line for all 


three thermal ranges are shown in the following Figures 6.2 to 6.4. 


C. "CONCEUSTONS 

From the calibration curves given in the previous chapter we can see that the 
portion of the curves in the temperature region used for the calibration procedure ts a 
straight line. The linearity of that region was the confirmation that the lrnear relation 
we Were getting was reasonable. 

The corrections we were getting from these straight lines were accurate to within 
an error of about 5%. The thermal values measured br the thermovision were 
accurate to within half a thermal unit. This caused an error of 1.5 % in the calculated 
temperatures. The variance of the best fitting starght lme for each thermal range caused 
a 3% error tn the thermal range 2, 3:5 % im the thermmdlirancess and 3. 7@7aemeune 
thermal range 10. Therefore, the total correction error was within 4.5 % in the thermal 
range 2, 5 % in the thermal range 5 and 5.2 % in the thermal range 10. Of course we 
have to keep tn nund that these caltbation functions are only valid for the temperature 


region about our target temperatures ( 9 to 18°C ). 
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‘ 4) 
CALIBRATION TEMPERATURE MEASUREMENTS 


TEMPERATURES IN DEGREES C TEMPERATURES IN DEGREES C 
BLACKBODY__ AGA DIFFERENCE BLACKBODY__ AGA__ DIFFERENCE 
EXACT COMPUTED EXACT COMPUTED 
3.0 16.1 13.1 15.3 26.9 i 
12 16.5 ee 16.0 27.3 - 11.3 
5.0 17.2 12.2 16.5 27.6 [1.1 
5.6 17.8 12.2 ene 28.0 11.0 
6.3 18.2 11.9 ig 28.3 10.8 
7.0 18.7 11.7 18.0 28.9 10.9 
8.0 19.7 ity 18.5 29.4 10.9 
8.5 20.1 11.6 19.0 29.8 10.8 
9.0 20.5 1.3 19.5 30.2 10.7 
9.5 21.2 11.6 20.0 30.7 19.7 
10.0 21.8 11.8 20.5 31.1 10.6 
11.0 22.5 1.5 21.0 31.6 10.6 
12.0 oer ties 21.5 32.0 10.5 
12.5 23.8 11.3 22.0 32.4 10.4 
13.0 24.5 11.3 os 33.0 0.5 
13.5 25.1 11.6 23.0 33.6 10.6 
14.0 25.7 ie 23.5 33.9 10.4 
l4.5 26.1 11.6 
MEAN DIFFERENCE _ : 11.03 
STANDARD DEVIATION : 1.724 
VARIANCE MSS 


Petes MRAIGHT LINE: Y = (1.0258) 4 = 11-62 
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TABLE 4 


THERMAL RANGE 5 
eee Oy Te viene URE Vieasl REMENTS 


Re wirehe he RES 1X DEGREES C TE VIPERA RES PD Girne eS 7c 


BLACKBODY Pome tir RE NCE BLACKBODY ANGE, UIPEEWeENCE 


Boe le CONPL TED EXACT COME TED 
0.0 14.1 14.1 16.0 Zo 10.7 
1.0 oI 14.1 17.0 ye 10.6 
ja 0) foro ee 18.0 Dome 10.2 
3.0 16.8 Ean 19.0 2o86. 10.6 
4.0 6 oe 20.0 0hs 10.5 
o.0 18.2 [oe2 21.0 31.0 10.0 
6.0 18.9 ie? EO neo ee 
7.0 he: 2 een) Be.0 9:6 
8.0 20.4 12.4 24.0 34.2 12 
Al oie 27 20) 34.8 9.8 
10.0 23 lS 26.0 cow ie 
bee0 23.0 Hee) Ze) Boo oe 
[20 24.2 12 28.0 31/26 9.6 
13.0 0 12.0 Zo) o0eo Ps 
14.0 25.7 ei 


mIEAS DIFFERENCE oS 
Se weRO DEVIATION. < o1/551 
VARIANCE gOS 


Peot pawns) RAG] GENE: Y = (1.1956) x - 16.6346 
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THERMAL RANGE 10 
Gye wie~ PEVPERATURE MEASUREMETNS 
eee EC RES Re DEGREES C eevir eee ese pierce eS 


BLACKBODY Pore UIFFERENCE BLACKBODY AGA JDIPrE RENCE 


Exel) COMPUTED EMAC] AGOuIPCIED 
0.5 14.2 (ee 16.3 27.2 10.7 
1.5 15.0 13.3 17.5 27.5 10.0 
0G 15.6 13.1 19.5 29.5 10.0 
aS 16.8 13.3 20.5 30.6 10.1 
4.5 7 13.2 DS 31.2 9.7 
5.5 17.9 12.4 aos 32.4 9.9 
6.5 18.6 12.1 23.5 32.8 9.3 
7.5 19.4 11.9 24.5 34.3 9.8 
3.5 20.1 11.6 p56 35.0 9.5 
9.5 20.6 Hl 26.5 aa 9.1 
10.5 21.0 10.3 27.5 36.9 9.4 
11.5 22.9 LL 28.5 mS 9.0 
12.3 24.1 11.6 29.5 38.6 9.1 
13.3 25.2 11.7 
14.5 25.8 11.3 

NES WIEFERENCE — : 11.00 

STANDARD DEVIATION : 1.489 

VARIANCE no 


PeseerltieslRAIGHT LINE: Y = (1.1807) X.- 15.656 
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VII. SIGNATURE MEASUREMENTS 


A. GENERAL 
All signature measurements were taken during mid May 1986, on the following 
dates and time periods. 
1. Mav 14 09:00 - 09:25 


Peewee 1S 11:37 - 13:07 
Cue eee 005 - 00:13 
peeeevidy 2) 17:18 - 17:37 
5. Mav 22 00:48 - 01:03 


Thermal imaging data were collected for the spectral range 8 - 14 fim using the 
LW detector of the AGA Thermovision 780. The thermovision system and _ its 
supporting computer Were set up in a laboratory building at the eastern side of the 
HOPKINS MARINE STATION grounds, at a window facing directly over the shore 
line toward the closest buoy. The laboratorv is about 8 - 10 meters above sea level. 

The target was the oceanographic ship R/V POINT SUR which was on an 
operational oceanography student cruise. During those periods in which it was made 
available for signature measurements the ship maneuvered in the vicinity of the off 
shore ( 635 meters ) buoy, displaying various aspects as requested. The distance 
between the the ship and our detector was about 650 m. At that distance, the closest 
the ship could reach, using a 79X70 FOV lens we had an image covering almost 2'3 of 
lenscrcel solwatmermocouples Installed about the superstructure were used to have a 
calibration of the actual ship temperature. The thermocouples were sampled at 
10-seconds intervals and the output voltages were recorded, to be later converted into 
temperatures. 

The environmental conditions for the time period in which the measurements 
Were takelare shown in the data presentation tables. Generally for these periods the 
area was under the influence of a high pressure system, with North West winds and 
coastal fog, especially during the morning hours. 

Thermal images for different ship aspects, were recorded and stored during time 
periods in which the ship was available for us and the environmental conditions 


allowed. All data were recorded assuming the ideal conditions appropriate to the 


Thermovision system measurements, that the target was a blackbody source and 
atmospheric attenuation was neglected. 


Fig 7.1 is a map of the Monterey Peninsula showing the area where the 
signature measurements were taken. 
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Figure 7.1 Area where the signature measurements were taken. 
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B. PARAMETERS FOR TEMPERATURE EVALUATION 


As mentioned in Chapter five the processing of the recorded and stored thermal 


images was done using the DISCO 3.0 computer program. For temperature evaluation 
the program asks for certain parameters to be provided. These parameters. as we have 
discussed before in the derivation of the exact measuremient formula. in Chapter five. 
are the atmospheric (air) temperature, the ambient temperature. the atmospheric 
transnuttance and the target emissivity. 

1. Atmospheric Transmittance 

The atmospheric transnuttance was calculated using the LOWTRAN 6 code. 
For our experimental conditions. temperatures about 10 to 15°C, we selected the 
midlatitude winter ( 45°, January ) atmospheric model for a horizontal path. The Navy 
Nlaritime was selected for the boundary laver aerosol model, discussed in detail in 
chapter four. For the geometrical path we inserted target altitude 10 m and-horizontal 
path length 630 m. The spectral range was set at 8 - 14 fim ( 710 - 1250 cm’! ) with 
frequency increment 5 en! 

The meteorological data needed, mainly wind speed and wind direction, for 
determining the air mass character, continental influence, were taken from 10-minute 
logged meteorological data recorded from instruments stationed on the ship. Such data 
eben Wrccotced for the last two observation periods on May 21 ( 17:18 - 17:37 ) and 
icwe #2 OO:48 - 01:03). For these periods meteorological conditions were 
approximated using observations taken from the meteorological station at POINT 
Pimms, tie closest station in the area of the experiment ( see Fig 7.1 ). Since the 
direction of the wind was from North West the continental influence as we can see 
from the map in Fig 7.1 was not very strong. So on a scale from 0 ( open ocean ) to 10 
( maximuni continental influence Vavie (senied mass caaincter Scale drom 2 10 57 

In Table 6 are shown the meteorological data used, the air mass character 
Seqiem wlG@sll ) and the calculated atmospheric trannuttance for each observation 
period. 

eepSSIV ICY. 

The emissivity of the ship surface was found by evaluating the emissivity of 
two metallic pieces taken from different parts on the ship. These specimens were 
painted with slightly different white paint, representing the paint coatings of the ship 
surface. These had also been exposed under the same environmental conditions for the 


Same time period as the whole ship surface. 
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INSERTED MIETEGR@2© 
TRANS viii NCE 
DATE VeIND SPEED? Aaihae eee 

MAY 14 (09:00-09:25) 4.80 ms 290 - 300 
MAY 18 (11:57-13:07) 6. 2051s 27/0 - 280 
MAY 21 (00:05-00:13) 5.20 mis 2/0 a0 
MAY 21 (17:18-17:37) 10.20 mis 520 

MAY 22 (00:48-01:03) 7 as 290 - 300 


DARE SG 


GICAL, DA TAsAN DA] VIOSPHEEG 


ICSTL ATMOSP? TRANSSE 


5 0.8677 
5 0.8539 
5 0.$640 
2 O.SL7S 
a 0.8282 


Assuming that the enussivitv was constant over the spectral region involved. 


we used Equation 5.5 derived in Chapter five, for our calculations 


= (1 a! 


Ty(lg - Ty) 


atm 


Making some more assumptions for the laboratory environment in which we were 


working, this equation was simplified further. 


First we had to assume that the ambient 


temperature and the atmospheric temperature Were the same in our case. This meant 


that the equivalent source thermal value for the surroundings ( I, ) had to be the 


same as the atmospheric path ( I 


O- latm 


Tea) 


atm. 


). So the above equation was then reduced to 


Second, for the distance of about 2 m between the detector and specimens being 


measured, we assumed no atmospheric attenuation, so T 


equation was reduced again to 
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& 


1, and the emussivity 


Where I’, is the Thermovision 780 measured thermal! value, L, 1s the calibrated 
thermal value for the known ambient temperature and [, 1s the calibrated thermal 
value for the known object temperature. 

With this last equation we could calculate the emissivity of the specimens. 
accurately enough. if we knew exactly the ambient temperature and the temperature of 
the specimens. [he procedure was started by heating the specimens as uniformly as we 
could. using a blower, to a temperature giving good thermal images above the ambient 
teniperature. This temperature was checked using a thermocouple. With the best 
thermal image on the screen. the temperature of the specimen, the ambient teniperature 
and the Thermovision thermal value measurement were recorded. The ambient 
teniperature was measured using a conimon thermometer. 

Knowing the exact temperatures, the calibration curves with an appropriate 
correction, discussed in Chapter six , were used to get the thermal values ( 1, ) for the 
specimen and ( I, ) for the surroundings. Therefore, we had what we wanted for our 
calculation. 

The above procedure was repeated several times at different teniperatures for 
each specimen. The results were averaged to the following values. 

a. 0.943 first specinien 
Gan U-205 Second Specimen 

These values are very close to the values we found in different books for the 
emussivity of all kinds of paint, which was 0.95. For the calculations done in this work 
We considered the emissivity of the ship surface to be 0.95, taking the average of the 
two values we found. 

3. Atmospheric Air Temperature 

The atmospheric air temperature measurements were taken from logged 
meteorological data recorded on the ship. The air mass temperatures measured close to 
our station had strong continental influence, so that they were higher than the 
temperatures measured at the ship location and not representative for our situation. 

The atmospheric air temperatures are shown for each date in the data 
presentation tables. 

4. Ambient Temperature 

These temperatures were measured using a common mercury thermometer 

near the open window. where our scanner was operating. Their vatue was between the 


higher room temperature and the lower outside temperature. 
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The ambient temperature fon eaeh date 1s Shown im thevG@diaeereseatarien 


tables. 


C.. TEMPERATURE DATA EVALUATION AND Teese NTA TiO. 

Inserting the correct parameters, the computer subprogram Evaluate vnage was 
used to evaluate the temperature distnbutiom Of the ships ior each inermaielemme 
recalled from storage. This temperature distribution was presented on the computer 
screen with the standard 8 color scaling combination. 

For better presentation of these temperatures we divided the ship image for each 
aspect in small squares and we tried to assign an average temperature to each square. 
‘AA cartesian svstem was used to present these temperatures for the starboard side, the 
port side and the stern. These temperatures were averaged for each aspect, for each 
date, in order to have the final temperature distribution for the particular aspect at the 
time period our data were averaged. 

The problem was that the aspect shown by the ship was not always the same, 
because she was drifting and had to keep in station using the engines. So pictures, say 
of the starboard side, taken under several time periods had different aspect angles with 
the line of sight. This caused reflection problems, with thermal radiation from hot 
surfaces such as the stack reflecting from cooler surfaces. Therefore, our scanner 
measured not only the emitted thermal radiation but both the reflected and the emitted 
thermal radiation from that surface. So the temperature measurement we were 
evaluating was higher than the correct temperature corresponding to that surface. To 
nuninuze this problem we tried to average the temperatures only from those pictures 
showing the most perpendicular aspect to the line of sight. 

However the temperaiures obtained were not correct, as discussed in detail in 
Chapter six, since the computer-calculated temperatures had still to be corrected using 
the empirical relation already derived. Hence the average temperature for each square 
had to be corrected using the best fit straight hne equation corresponding to the 
thermal range in which the data were collected. 

The corrected average temperature distribution for the starboard side, the port 


side and the stern, available for each date, are given in the following Tables 7 to 11. 
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D. TEMPERATURE DATA ACCURACY 
As mentioned before, six thermocouples were installed about the superstucture, 
at the following locations: 
- Tl: On the after starboard side of the deekhouse 
- 72: On the forward starboard side of the deckhouse 
-  T3: Qn the aft bulkhead of the bridge ( away from the stack ) 
- T4: On the aft bulkhead of the deckhouse 
- 75: On the port side cf the deckiouse 
- [T6: On the after deckhouse forward of the stack 

These thermocouples were sampled at 10-seconds intervals and the output 
voltages were recorded and later ¢corvetted Wito tenioeri nes aus, We Shame 
indication of the actual temperature at six different locations on the ship. for the whole 
time period we were taking our measurements. In this way we had an absolute 
reference to conipare our data with to estimaie the “aceuracy of (thes pioces cn 
calculating the temperature of the target from measurements of thermal radiation. 

First we had to ht the thermiocouples im the squares covering the areas aicn 
they were located. Then we had to find the mean calculated temperature corresponding 
to these locations. For example for a thermocouple located at the intersection of two 
or four squares, the calculated temperature corresponding to this thermocouple was 
taken as the mean temperature of the two or four squares. Finally we had to compare 
the above temperature with the thermocouple temperature averaged over for the same 
period as the calculated temperature for each square. 

In the following Table 12 are shown the thermocouple locations and the 
comparison of the actual with the measured and calculated temperatures in degrees C 
for each date. 

As we can see, for most of the time the temperature dillerenees were sitll Pie 
big negative differences in the stem drea for the twoniist dares (Siemperauineses 
thermocouple higher than AGA thermovision temperature ), are believed to occur 
because the thermocouples were not yet stabilized at that period. That 1s, at a time 
when the ship was turning, hot air from the stack area was blown onto these 
thermocouples, making their measurements fluctuate. So for measurements at this time, 


just after the turn, the thermocouples did not have the time to be stabilized. 
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TABLE 12 
THERMOCOUPLE - AGA THERMOVISION TEMPERATURE COMPARISON 





w STARBOARD @® PORT + STERN 


Ate MAY 14 (09:00 - 09:25) fA 28 (22:57 - 13:07) 

et £0 ate ed 

SENSOR S==sMoc. AGA a Se LeeeRMOC. AGA DIETER. 
fea) = Boe 3 a. * A .- a =a ~ 
~ = > _ e = 7 38 = wd « 35 me 84 -_-_ $2 a ¢G 08 
x 
= oe! g2. 42 EZiuaoe = 22 L226 12.84 = <2 
by ~ en ed Aa 7 A = 2 = = 
w ao 25. 4. 22.23 - 1.26 eee hs —Z. 65 - 1.47 
oO 
Ry 
Fe Le $3.27 33.02 = 6.25 14. 902 ee eyes = 4.22 
fz] rf - - = ~ * « * a «= 
ce —) Lay, 14.07 =2:2: 20 Loe oe ize! a2 630 
Y 

t 15.84 eS - 1.65 14.95 LZ 80 = 92.15 
DATE= MAY 22(00:05 - 00:13) 
Sanson ESERMCC. AGA we So. 


=} 12.06 22. 55 - *.53 


eZ Loess 2 ah S Soo 2 


STARB 


YATE Mey 23 (17:16 - 17:37) Mow 22°5(00:45 = 02:03} 
Ser5OR THERMOC AGA DIFFER THERMOC AGA DIFFER. 
oe we.28 14530 + 1.02 52 9,99 - 0.47 
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E. RADIANCE CALCULAIIG. 
Knowing the temperature distribution of the ship, we could calculate the radiance 
for the spectral region of the measurements. 


Starting With the Planck blackbody radiation law we could write 











h Is 
Lt | oanhe? W 
M = J W(A,T)dA = | (———) (————_) dd 
> ete ee 
i.» h5 d on Ol | cm 
5 he 
Lf we Let ae eee cr ea we could write the above equation as 
follows K 
hy Cy W 
M = {| —-—_---_—_ dh 
we Ne (eC AT - 1) em? 


Solving this integral in the spectral region 8 - I4 ttm. for the Known temperatures 
we found the in-band radiant flux density M in W cm*. With the assumption that our 
Surface was an opaque Lambertian surface, we could divide the radiant flux by m and 


multiply by the emissivity ( € ) to obtain the required spectral radiance, 


£ VI W 





Z 
1 Cilesl 


The integration was done numerically using a program for the Tl - 59 calculator. 
details of which are given in Appendix B. The results of our calculations are shown in 


the following Tables 13 to 17. 
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Vif. DATA ANALYSIS 


A. INTRODUCTION 

In this chapter we will analyze how the sea background radiance was affected due 
to reflected sky thermal radiation. We will “discussm@iow sthe reflectance “nde 
emissivit¥ of a sea surface roughened by various wind speeds atfected the background 
radiance. Also we will calculate the apparent temperature of the sea surface, and will 
find the effective target to backsround temperate wdimerence sensed Dy tire 


Thermovision under the specific environmental conditions experienced. 


B. EFFECTIVE SEA SURFACE TEMPERATURE 
1. Sea Surface Radiance 

As mentioned in Chapter III the infrared radiance of the sea is the sum of the 
radiance due to thermal emission and the radiance due to reflected incident radiation. 

The reflected incident radiation in the spectral region of 8 - I iin 1ssonim 
thermal sky radiation bv dav and might. The contribution of this radiation to the total 
sea radiance depends on the reflectance of the sea surface (R). 

Because the water 1s éssentially opaque ineine above spectralereorone ric 
reflectance also controls the emissivitv of the sea surface, since for an opaque matenial 
we can writeeé = 1-R. 


Therefore, we can-wtite tat the total ridianes Giitic sea’ sUriices. 


Lo=(1-R)L 


W 


a obs Be : (eqn 8.1) 


Where L, 1s the radiance of the sea surface. L,. is the radiance due to sea water thermal 
emission and L,,,. 1s the radiance of the sky. To write this equation we assumed that 
the emissivity of the sky, as we have previously discussed in detail, is one. 
2. Reflectance of the Sea Surface ; 

An important factor in determining the reflectance of a surface is its geometry. 
For the sea surface this geometry depends on the slope distribution of the waves at 
various wind speeds. 

In our discussion about the geometry of the sea surface we have shown, in Fig 
3.8, the difference in reflectance between a flat sea surface and a surface roughened by 


a Beaufort 4 wind speed. 


These results, as Fig 3.8, were obtained by Cox and Munk, working in the 
visible range with solar radiation. Such information is tacking for the infrared and night 
conditions, although the above curves seem to be helpful in this area also. We make 
this assumption since the reflectance of the sea water at zero angle of incidence is 
almost equal to the average reflectance in the infrared region from 2 - 15 ttm ( see Fig 
3.4 ), and the sky radiance variation is alniost the same for dav and night conditions. 

For our analysis we used the same curves corrected for the average index of 
refraction in the spectral region $ - 14 ttm. Actually we wanted to find a correction 
factor for the reflectance. Since we knew the index of refraction of water for both the 
visible and 8 - 14 jim regions, we could find the corresponding reflectances at normal 
incidence using the equation R(O) = (n- 1 ea easel >. Their ratio gave us the 
required correction factor. This factor was used to correct the given curves. assuning 
that the correction we found for normal incidence applied for all incidence angles. 

[Net 7)Seisies the index of refraction used to caleulate the sreflectance. for 
uncontanunated sea water, as n = 1.338 and quotes the corresponding calculated 
value of R,(0) as 0.020. Even though the reflectance equation for n = 1.338 gives 
R,{0) = (0.02089, we will use for our calculations the valueof 0.020 used to draw the 
curves. The average index of refraction for the region 8 - 14 ttm was found from Fig 
Peet omec ie —) 1503 > with this index of refraction the reflectance was calculated as 
R; (9) = 0.01731. As we can see the reflectance for the region § - 14 ftm at normal 
incidence is lower than the reflectance in the visible region, in the ratio of R:(0) R.(0) 
= 0.865. With this correction the curves were shifted a little lower. 

The following Tables 19 to 23 at the end of this chapter show the reflectance 
ie valeseame tne enussiity ( 1 - R’) values as they were found from) the corrected 
curves, for a flat sea surface and one roughened by a Beaufort 4 wind, for certain 
angles of incidence. Here we reiterate that the values obtained for mcidence angles near 
the horizon are less accurate, because the curves are not strictly applicable in this area, 
aomeascusscammecmapter 111. 

3. Apparent Temperature Calculation 

As we said before in section one, the sea surface radiance can be written as 
taewOllowine sum ( Eqns.! ) 

Peg Ri) Laat R Lexy 


that 1s, that the radiant flux of the sea surface is the sum of the emitted radiant flux 
and the reflected incident radiant flux. For the whole spectral region we can write 
oe 4 4 > 
of. (1 oa aa ROT 44, ; (eqn 8.2) 


where T, is the apparent sea surface temperature, T,, is the actual sea surface 
temperature and T,,.,. is the temperature of the atmospheric air. Dividing both sides of 
the equation by 6 we get the more simphfied form 

T, =(1-R) 1,7 = R Tyee (eqn 8.3) . 

With this equation we could calculate the apparent sea surface temperature, 
the temperature that our thermal svstem would sense if it were operating in the whole 
spectrum region and the atmospheric attenuation was neglected. 

With the above assumption, we calculated with moderate precision. for the 
region of observation, the apparent sea surface temperature for a flat sea surface and 
one roughened by Beaufort 4 wind. Calculations were done for incident angles at 89°. 
where we had been working, 85°, 80°, 75° and 70°. The reflectance (R) for the two 
kinds of surfaces and the incidence angles listed was found using the curves in Fig 3.8. 
after being corrected for the spectral region 8 - I4jftm. The actual sea and air 
temperatures were taken from the recorded meteorological data on the ship. 

The results of these calculations are shown in the lables 19m6 25 at theveaa 
of the chapter. Detailed analvsis and discussion of these results is given in the following 
section. 

4. Calculated Temperature Data Discussion and Presentation 

From the equation we used to calculate the apparent sea surface temperature 
We can see that the two terms at the rlgmt side are 1 ercerm Telated Tomme wercetalie mem 
the sea surface. 

Therefore, when the reflectanee of the sea suriace is hich, the term Wwitneeue 
sky temperature ( T.,,, ) dominates, and conversely when the reflectance is low the 
term with the actual sea surface temperature ( T,,. ) dominates. So whenever the 
reflectance was high the apparent sea surface tenyperature we had calculated was closer 
to the skv temperature. Since for our case the sky temperature was always lower than 
the actual sea surface temperature, the apparent sea surface temperatures we have 
calculated were lower than the actual. Om the contrarywiicie cme oleeraniece: as 
low, which means that the enussivity of the sea surface was high, the apparent sea 


surface teniperatures we have calculated were elese siominemacmalmescamestini aes 


$6 


temperatures, but still some tenths of a degree lower, because of the smiall contribution 
from the second term. 

Looking at the reflectance values given for the two kinds of surfaces, at 
Various incident angles, in Tables 19 to 23 we see that. for large incidence angles (i.e.. 
grading incidence) the difference in reflectance of the two surfaces is large. But as the 
incidence angle decreases the difference decreases and at about 70° goes to zero. (for 
miore details see Fig 3.8). For this reason the difference between the apparent 
temperatures we have calculated is large at high incidence angles ( 89°, 85° ) but goes 
to a few tenths at 70 °. 

The calculated temperatures are shown in the T, row of the Tables 19 to 23. 
In the Ty row the difference between the actual sea surface water temperature and the 
calculated sea surface temperature is shown. T,, 1s the temperature computed by the 
Thermovision 780 system, divided by the atmospheric attenuation. This was done 
because we wanted to compare this temperature with the T, temperatures. But since 
the IT, temperatures do not include any atmospheric attenuation, it was necessary 
either to include the atmiospheric attenuation into these values or amplify the 
Thermovision measurement, to obtain a value without the atmospheric attenuation. 
We followed the second method and we got the values T,,, shown in the Tables 19 to 
2), 

Comparing these T,, values with the T, values we have calculated for incident 
angle at 89°, we can see clearly how much the reflectance of the sea surface, 
roughened at various wind speeds, was affecting the temperature sensed by the 
Thermovision. 

Finally in the T,,4 row Oleic: Wael ssl tome wetinc m@itercnce between, the 
actual sea surface temperature and the T_,, temperature is shown. As we have expected 
this difference is big at low wind speeds when the reflectance 1s high and the apparent 
fem@@merature 1s Closer tO the sky temperature. 

Generally we can say-that the apparent sea surface temperatures found were 
lower than the actual sea surface temperatures, since the skv was always colder. Their 
difference was related to the incidence angle and the roughness of the sea surface, with 


higher values at incidence angles near the horizon and lower wind speeds. 
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C. TARGET-TO-BACKGROUND TEMPERATURE DIFFERENCE 
1. Effective Temperature Difference 

An important factor for the target detection probability of an infrared svstem. 
is the target-to-background temiperature difference. The higher the temperature 
difference between a target and its background the greater the probabilitv and the 
distance of detection. 

But as we said before the effective temimeratitire diierence ( AT) taco 
thermal systen. would sense, under various environmiental conditions. looking at 
various incidence angles, might not be the samie as the actual temperature difference. 

Therefore, the precalculated detection distance, using the actual temperatures. 
at which We expect our system to detect a specific target, may not coinside with the 
actual distance at which the target is detected by the system. 

In this analvsis we attempted to find the effective AT, sensed by the 
Thermovision, under the various environmental conditions worked in, and compare it 
with the actual AT . The effective AT was found by taking the dillerence benveenmtie 
temperatures the system calculated for the hottest spot on the ship, ie. the stackwac. 
the sea surface. These temperatures, were calculated under the ideal conditions in 
which a thermal imaging system operates. considering emissivity one ( € = 1 ) and 
neglecting atmospheric attenuation ( t, = I ), and.were corrected using the calibration 
relation we had derived. The actual AT was found using the temperature of the stack. 
computed from the radiation measurements using the exact parameters (emissivity, 
atmospheric attenuation, etc. ). 

These teniperatures, after being corrected for cahbration purposes. were the 
most accurate available for our case, since we did not have any thermocouples on the 
stack to give a direct indication of the actual temperature. The temperature of the sea 
surface was taken from the meteorological data recorded on the ship. 

From this analysis we found that the effective AT, for the environmental 
conditions of the measurements and incidence angles near the horizon. was always 
higher. by about a factor of 1.6 to 2 .than the actual AT. This was an indication that 
Variations in enussivity and reflectance of the sea surface under the experienced 
environmental conditions had affected the effective AT, sensed by the Thermovision. 

Our results are shown in Tables 19 10°23 ar thesend@ ol the chaprer: = mene 


detailed discussion of these results is given in the following section. 
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2. Temperature Difference Data Discussion and Presentation 

The AGA temperatures shown in the Tables 19 to 23, for the hottest spot on 
the ship and the sea surface are temperatures calculated by the computer program 
using the thermal values measured with the Thermovision 780, assuming € = | and t 
= | and corrected for calibration. 

Therefore, these temperatures correspond to the thermal radiation attenuated 
by the atinospheric effects and the emissivity of the measured surface. Thev are the 
effective teniperatures our svstenl was sensing under the environmental conditions 
experienced on each date, at high incidence angles near the horizon ( about 89° }. 

The values we present in the Tables 19 to 23 for both the effective and actual 
temperatures are the average values from the data we obtained for each observation 
date. 

Comparing the effective with the actual AT values we got for each observation 
period, we see that the effective AT for the environmental conditions experienced was 
always higher than the actual. It is obvious that there are two main factors causing the 
MioMenic eet rs |, tlhe emissivity difference between the two surfaces and the reflected 
wicident radiation. Since the enussivitv difference was small. the dominant factor was 
Pveereliccted thcident radiation. Because the reflectance of the ship surface was too 
pmiatt(e — 0:95) tomcause such a difference. it was the reflectance of the sea surface we 
had to check. But we have already done this in the previous sections the calculation of 
the apparent sea surface temperature. Then the conclusion was reached that the 
apparent sea surface temperature was lower than the actual sea surface temperature 
because the sky was alwavs colder than the sea for the local conditions. These lower 
teinperatures Were related to the inctrdence angle and the roughness of the sea surface, 
with the lower values noticed at incident angles near the horizon and low wind speeds. 

Hence the Thermovision was sensing lower sea surface temperatures due to 
reflected thermal radiation from the colder sky. This was causing the effective AT to be 
higher than the actual. Therefore. we would expect these lower sea surface 
temperatures and high ATs to be related to the wind speed conditions which control 
the roughness and consequently the reflectance of the sea surface. stnce the incidence 
angle was effectively constant (near the horizon. about 89” ). 

From the tabulated and plotted (Fig 8.1 at the end of the chapter) ratio of the 
effective and the actual AT for the various wind speed conditions we experrenced. a 


significant trend can be seen. Even with the assumptions and cumulative errors 


S9 


inherent in the approach it can be seen that for the lower wind speed of 4.8 mson the 
14th of May we had an effective AT higher than the actual by about a factor of 2.01 
and for the higher wind speed 10.2 m/s on the 21S of May ( 17:18 - 17:37 ) we had an 
effective AT higher only by about a factor 1.63. 

Thus as expected for the lower wind speed) higher reflectance, cunigeum: 
higher AT and for the higher wind speed, lower reflectance, lower AT. It was obvious 
that the effective AT was strongly dependent on the reflected incident sky radiation. 
Also within the spread of the plotted points we may conclude that the AT ratios are 


proportional to the wind speeds. 


D. CALCULATION EXAMPLE 
1. Introduction 
To help understanding of the following Taslesm1 0 10) 23a) emer cure ccue 
catculation example to show how the data presented were obtained. In this example we 


h 


use the information we had for the 18"! of May experiment. 


DATE-TIME: 18_@MAY  11:57- 13:07 


PMY DeMidelweNl is (on) oa 
SEA SURFACE Te ytV ERE iitiNE (1) meine 
WIND SPEED - DIRECTION : 6.26 m's(B-4) 280° 


OVERCAST, HAZ), awe | ye New 
ATMOSPHERIC TRANS MILDANG@E nee? 


Average temperatures computed from Thermovision values for the hottest 
point on the ship (Hot Spot) and the sea surface temperature (S.S.T.), (assuming ¢€ = 


land t, = 1) are as follows: 
UNCALIBRATED VALUES : HotSpot : 22.80°C 
SSslaee GG 


CALIBRATED VALUES > Homscom lie 2aG 
Sle oe © 


Alt data for this date were collected using THERMAL RANGE 2 so the calibrated 


values were found using the straight line equation Y = (1.0238) X - 11.62 . 


90 


The relation between the Beaufort sea state scale and wind speed (m,s) is 


shown below in Table IS. 


WepLeE ls 
Dee | OR] SCALE 


B WIND SPEED enaes 
0 0.00 - 0.46 
| O.51- 1.S4 
2 2.06 - 3.08 
3 3.60 - 5.14 
4 SGD co. eee 
5 §.74 - 10.80 
6 Pi = 1aase 


2. Apparent Sea Surface Temperature 
To calculate the sea surface apparent temperature we used Equation 8.3 
4 . 4 4 
alk = {( ome lee ste RT ky 


Using the values in Fig 3.8 for incidence angle 89°, where we have been 
working, we got the reflectance values for the flat sea surface and that roughened by a 
Beaufort 4 wind. These values, corrected for the spectral region § - 14 f’m_ using the 
correction factor R,., Ry. = 0.865, are the following: 

FLAT SURFACE R = 0.804 
ROUGHENED SURFACE” R=. it 


Therefore. from: Equation $.3 we get for the flat surface 
i OOO temas een Ost 109.7 2 
T, = 10.65°C , 


and for the rough sea surface at Beaufort + wind speed we get 
ie: = GeroniO Gleam e ea )atOn O17). 
eet eC, 
The Tg values represent the difference ( T,,. - T, ) between the actual sea 


surface temperature and the apparent temperature we calculated. [lence, we have for 


oi 


the flat surface 13.1 - 10.65 = 2.45°C and 13.1 - 12.84 = 0.26°C for the roughened 
surface. | 

The T,,, value is the sea surface temperature computed by the AGA system. 
but without the atmospheric attenuation. To get this value we proceeded as follows. 
The uncalibrated value which the svstem calculated for the sea surface, 20.65°C., is 


aniplified by the atmospheric transnussion T,, = 0.8539; hence, we have 


¢ 


20165 = 0.8539 = 2440°¢ 


This value was corrected for calibration purposes using the straight line equation for 


thermal raricceeato: ive 


Ty = 24-40 ( 1.0238 ) - 11.62 = 12.83°C. 


As we would expect, this value must be close to the T, value for the roughened surface 
( B-4 ), since the wind speed on this date was 6.26 m.s ( B-4 ). 


The eae values are the differences eel ean indicate how much colder 


in 
the effective temperature sensed By the Thermovision was than the actual sea suriace 


temperatures hlencemuc tiawe 


Tean= 1310123 sane ze C. 


5. Effective Sea Surface Temperature 
As we have said before, the values given for the AGA are the values the 
svstem calculated, corrected for calibration ~ Phese valtes are mot corrected) tomeca. 
atmospheric attenuation and emissivity. Therefore) for the hot spet we haveum. 


uncalibrated value 22.8°C corrected to 11.72°C, using 
AGA HOY SPOT = 22.8 (i¢ess ieee — 1172 C 


Simularly, for the value of the AGA sea surface temperature we have the uncalibrated 
value 20.65 °C which corrected gives the value 9.52°C. 

For the actual conditions, the temperature of the hottest spot on the ship, 
since We did not have any thermocouples on this location . was taken as the highest 
temperature we found for this location in the calculated temperature distribution of the 
ship. This value for the 18" of May was 14.28°C. The sea surface temperature was 


taken from the recorded data on the ship. 


oe 


The AT value represents the hot spot to sea surface temperature difference. 


Hence, we have effective AT, = 2.20 °C, from 
en see = 220 C 
ie actual a — elec. from 
AT, = 14.28 - 13.10 = 118°C. 
The ratio AT, AT, = 1.864 represents the fractional increase in effective temperature 


for the particular environmental conditions encounered. 
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IX. CONCLUSIONS AND RECODMMENDATIONS 


A. CONCLUSIONS 

Temperature and radiance distributions of an instrumented target measured 
remotely by radiometric imaging with the AGA Thermovision 780 system have been 
compared with locally measured environmental parameters. The ship-to-background 
temperature difference, frequently used as a measure of the contrast. was compared as 
measured thermometrically with local temperature sensors and radiometrically using 
the Thermovision. Significant differences have been observed, depending on the time of 
the day, windspeed, and atmospheric conditions. These differences have been assigned 
to effects of the low emissivity and high reflectivity of the sea surface at near grazing 
angles. 

Reflected thermal radiation from the sky contributed greatly to the sea surface 
radiance sensed by the system. Comparison of target-to-background tempcrature 
differences from radiance and from thermometric measurements showed large 
discrepanctes due to the effect of the reflected sky rathation, Characteristic Ol mine 
atmospheric temperature. This reflected radiation was strongly dependent on the 
roughness of the sea surface under conditions of varving wind speed. The differences in 
the air and sea surface temperatures caused the effective (radiometric) temperature 
differences between the target and background to deviate both positively and negatively 
from the true temperature difference. With the skv colder than the sea surface, as was 
always found in the first series of measurements, the reflected sky radiation caused the 
sca surface temperature measured by the Thermovision to be lower than the actual. As 
a result, the effective temperature difference as observed with the Thermovision was 
higher than the actual difference. Both the temperature decrease and the increase in the 
effective temperature difference were related inversely to the sea state conditions. 

An expression for the sea surface radiance including the reflected sky radiation 
has been developed and used to model the radiometric effective ship-to-background 
temperature difference. This model can be evaluated for known atmospheric conditions 
using the LOWTRAN 6 computer code and has been found to be consistent with 


observations to date. 


100 


In analysis of radiometric images of the ship taken at different aspect angles, 
apparent differences in radiance distribution were observed. This is considered to result 
from thermal radiation from hotter portions of the ship reflecting from cooler regions. 
increasing the apparent surface temperature. Consequently, temperature values from 
only those surfaces more nearly perpendicular to the line of sight were averaged to 
obtain representative temperature and radiance distributions. 

Although the quantitative results cannot be considered to be of great accuracy. 
they show clearly the strong dependence of the reflected sky thermal radiation on the 


Sea roughness due to the wind speed at the times of the observations. 


B. RECOMMENDATIONS 

Since we had to calibrate our system empirically to achieve the accuracy of our 
mieasurements it is Obvious that we have inserted some error factors into our results. 
Therefore. we recommend a new series of target and background measurements to be 
taken with the whole system being accurately calibrated. 

A few months before this work was completed the Thermovision was tested and 
calibrated accurately. The gain of the preamplifier in the scanner was adjusted and new 
Silmratien cUrwes were derived. However there remained a problem in that the 
coniputer did not calculate the exact temperatures because the thermal values mserted 
by the svstem into the computer prograni were off bv four thermal units, even though 
a few weeks before this work was completed we took a new series of measurements of 
the same target. Our results and analysis are shown in Appendix A. 

All our measurements were taken at an incidence angle near the horizon where 
Paercurves 1) Fie 3.Squised togiind thesrellectange,are not strictly applicable. Therefore, 
depending on the operational use, we may need this information with greater accuracy, 
we recommend: 

For systems like missiles or airborne forward looking infrared devices a series of 
measurements should be taken for incidence angles between 70° and 80° where the 
curves in Fig 3.8 are applicable, and where more accurate results can be obtained. 

For forward looking devices operating on ships, with incidence angles near the 


horizon, long series of measurements, under all possible kinds of environmental 


conditions, may be needed to provide sufficrently accurate information. 
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APPENDIX A 
DATA ANALYSIS OF LATER MEASUREMENTS 


lL. GENER Ie 

A second set of thermal imaging data were collected in mid Novermnber,1986. All 
measurements were taken in the spectral region 8 -14 Itm as before. using the AGA 
Thermovision supplemented with a microcomputer. The ship target was again the R V 
"POINT SUR", scheduled on an operational oceanography Sitdentecruise, anda 
geographic location where the thermal imaging device and the ship were positioned 
during the experiment was also the same. 

These data were collected in order to obtain further, more accurate. 
measurements, under different environmental conditions, after the Thermovision had 
been calibrated. This was done to check our previous results and improve the 
information on the reflected sky thermal radiation and the effective target-to- 


background temperature difference relation. 


2, CALIBRATION TEST 

The Thermovision was accurately calibrated after the first measurements had 
been taken. The preamplifier gain of the scanning device which was found to be 
incorrectly set, was adjusted and new calibration constants were derived and inserted in 
the computer program. However. even if the Thermovision was operating accurately on 
this occasion. the temperatures calculated’ bv tie Cemputcr were not accumiie aime 
problem lav in the thermal values transmitted to the computer from the Thermovision. 
The thermal .waltes trccetjedmar sig computer were about 3 104 thermal winitssie. « 
than measured directly on the Thermovision. The error was minimized by adjusting 
the THERMAL LEVEL control to read about 4 thermal units higher, using the 
THERMAL LEVEL ADJ control on the Black Wy lite monitor chassis. In this wayweumte 
the computer received almost exactly the thermal values corresponding to the object 
in view. Of course this was not the night way to solve the problem, since now the 
Thermovision was measuring thermal values 4 thermal units higher than previously. 
However, the system was empirically calibrated again, following exactly the same 


procedure described im Clizpter . 17 
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The calibration temperature measurements for the two thermal ranges used in 
these series of measurements ( 2 and 5 ) are shown in the following Tables 25 and 26. 
The data were plotted as shown in the following Figures A.1 and A.2 and the best fit 
straight line was found. As can be seen from the mean and variance of the difference 
between the exact and calculated temperatures. the temperatures the computer 
eileWiitcdewere very close (oO the exact temperatures. The correction error found this 
tine Was Within 2.9 °> for the thermal range 2 and 3.5 °% for the therinal range 5. Of 
course again we have to keep in nund that these calibration functions are valid only for 


. . - 0 ~ 
temperatures in the region from 5°C to 25°C. 


3. MEASURED DATA 

Thermal images of the target and its background were collected on the jst) of 
November ( 09:07 - 09:12 ), om the 16" of November ( 17:58 - 18:02 ) and on the 18th 
of November (08:34 - 08:36 ). All measurements Were taken using the saine equipment 
as before, in the same spectral region and geographic location as is described in 
Sieipte@ewiles | te omy wilterence was that this time the R,YV POINT SUR*“did not 
maneuver near the buov to show different aspects but passed by the buov on her way 
into and out of the Monterey harbor. The distance between the lab station and the 
Ship was about 650 m on the 14¢) and 1st? when the ship was going out of Monterey 


Bav showing the port side and about 750 m on the por 


on her wav into the harbor 
showing the starboard side. 

The environmental conditions experienced during these measurements are shown 
ineeeieminiaepresemiatian tables for each observation period (Tables 27,28,29). 
Generally we had low speed winds blowing from the North West. overcast skies and 
fog on the last day. The air temperature was always higher on this occasion ( within 
0.5°C to 1.9°C ) than the sea surface teniperature. 

Wemanal.ze the daia the same procedure described in Chapter VIII was followed. 
The apparent sea surface temperature was calculated and the effective target-to- 
ClencrOUmemmicniperaiure difference, sensed by the Iherniovision under the local 
eo oumemn ceoncamions. was found. These résults are presented in lables 27,28,29 at 
the end of this appendix. Also a plot of the ratios of the effective to actual ATs versus 


Pimdespecd isioresented in Figure A.3 alter the tables. 


TABLE 


THERMAL RANGE 
NEW CALIBRATION TENIVEKAT 2 eee sig 


TESIPERATC RES ENE Gy Eicec 


BLACKBODY AGA DIE PEE Nee 
Ea | COMPUTED 

on) 4.7 0.3 

6.0 5.8 Or 

70 ore 0.3 

8.0 7.8 re 

220 Sa 0.3 

10.0 oa 0.3 

11.0 ie -Q. | 

Ie 12.2 -Q.2 

ae es 0.2 

14.0 Eee Ue 

15.0 14.9 0.1 

16.0 is 0.1 

£ae0 [fee -0.2 

18.0 18.2 -Q.2 

PoC is 0.1 

20450) 202 } 252 

0 29 Ox | 

BMAD Zon 0.0 

2500) 2320 0.0 

24.0 24.1 -Q.1 

2&0 Loe =iaz 
VIEAN DIT EERE \Ge : 0.0286 
STANDARD DEVIATION. 022505 
VARIANCE > 0.0506 


BEST FIT STRAIGHT LINE YY =(09793) <=] 0c2. 
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[AL o 
NEW CALIBRATION TENPERATU RE DATA 


TEXIPERA TOS eerie sc 


BLACKBODY AGA DIFFERENCE 
EXACT COMPUTED 
0) 4.6 O.4 
6.0 5.5 0.5 
7.0 6.6 0.4 
8.0 7.6 0.4 
9.0 or 0.6 
10.0 9.6 0.4 
11.0 10.5 0.5 
12.0 L1.d 0.6 
13.0 12.5 0.5 
14.0 13.6 0.4 
15.0 14.6 0.4 
16.0 15.7 0.3 
17.0 16.7 0.3 
Iga 17.7 0.3 
19.0 18.6 0.4 
20.0 19.6 0.4 
21.0 20.7 0.3 
22.0 Died 0.3 
23.0 22.6 0.4 
24.0 23.7 0.3 
25.0 2G 0.4 
VIEAN DIFFERENCE : 0.4047 
STANDARD DEVIATION : 0.0921 
VARIANCE : 0.0081 


BEST FIT STRAIGHT LINE ¥ = (0520) 
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4. DISCUSSION AND ANALYSIS OF THE RESULTS 

As was expected, reflected thermal radiation from the warmer skv affected the sea 
background radiance under the various sea state conditions due to the wind speeds 
encountered. Therefore, the apparent and the effective sea surface temperatures were 
found to be higher than the actual. As before. the differences between the apparent and 
actual sea surface temperatures were higher at incident angles near the horizontal and 
at lower wind speeds. 

Simee the effective sea surface temperature Sensed oy the [hermovisione a, 
higher than the actual, the effective AT Sensed bY the detector was lower tian 
actual AT. 

The major error factor occurred because we did not have an indication of the 
actual temperature of the stack, the hottest spot on the ship, and we had to trust again 
the teniperature the system calculated for this spot. On this exercise we did not have 
thermocouples positioned at various locations on the ship as we had before. to 
calibrate the measurements. An attempt was made to record a measure of the stack 
temperature, using a common thermometer fastened to the stack surface. 
Unfortunately since we could not place the thermometer on the surface viewed by the 
system, we did not get accurate measurements. 

Even with all the uncertainties about the exact ATs and their ratios, the results 
show again that the variations in enussivitv and reflectance of the sea surface under the 
locat environmental conditions had affected the effective AT sensed by the 
Thermovision. 

From the tabulated and plotted ratios (Fig A.3) of the effective and actual ATs it 
can be seen that for the lower wind speedmoi 1-5 aioe pth ~Nov.), and the higher 
reflectance, we find the lower AT ratio of 0.939 . This ratio increases to 0.954 when 


the wind speed increases to 2 11s { pth 


Nov. ) and the reflectance decreases. Normialiv 
one would expect to have the higher ratio on the 142 of November when the wind 
Speed was 2.9 ms, but unfortunately the measured data do not seem to be consistent 
for that date. The temperature value recorded for the stack on this dav appears 
inconsistent with the reliably measured values of air temperature. However, the ratio 1s 


close to one, mol too far from wnat Ws ex peereas 
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5. - CONCLUSIONS AND RECONIMENDATIONS 

It is noticable again in this series of measurements that reflected skv radiation 
affected the sea surface radiance sensed by our system. Since the sky was warmer than 
the sea surface , during the time periods of observation. reflected skv thermal radiation 
caused the sea surface temperatures sensed by the Thermovision to look higher than 
the actual sea surface temperature. The higher effective sea surface temperatures 
caused the effective ATs sensed bv the Thermovision to be lower than the actual ATs 
under the various wind speeds experienced. 

To improve the accuracy of the data set, it 1s recommended that a series of 
ednements ©e taken, swith the wiaole system correctl cahbrated and with a 


measurement of the actual temperature of the side of the stack viewed by the system. 
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TABLE 26 
SEA SURFACE PEMPERATURE AND EereC (Vee a 


DALE / Live ¢ ee NOV | OSs 7 ao ore 


© 


WIND SFEED 7 DIRECTION :) 2ygem/sec (o-2) eezo 


LIR TEMPERATURE : “Toor 

AMBIENT TEMPERATURE : 25°C 3 
SEA SURFACE TEMPERATURE : 12.2°C 

TARGET Diss aNes 650m 

OVERCAST, VISIBILITY 4 NM 

ATMOSPHERIC TRANSMITTAN ©. 8864 


| ANGLE OF INCIDENCE AND SUREAC= Conte cn 
| 

















| | | eS | es° | sc° | ce i 70 | 
| | | PLAT | p-4 |lenat | en) Gewese | leah SN ene Wee | 
IR lito. 804|0. 110| lo. 513| Cc. 238| | 0. 292\¢. 1665/0. 180) ¢ 5i| lo. 120|0. 120 | 
i-R |/0. 19616. sie 487|0. 86 2| 10. 708|0.834 |0. 820) 0. wale 880 | 0. 880 











| Tn = 225/56 TP 

E= | -0. 36 | 

| remper.| wot spot | sea surrace | At | RAT: | 
1 aca | Qewecmes | NNOReR | 3.87 Gemeanos | 4 
| ACTUAL | “Secor! || eer. | 3.70 | 


R > RESLESTANCE “O8 THE SEASURees 

ELAZ: EEA essheounence 

B-¢ : SEA SURPACERRCUGEIN=D (2 Ss205er eee 

Le ALCULATED eo Sia SUREAC =n oe oe ner 

fa : DIFSERENCE BETWEEN ACIUAL AND Galeuul bay oe eee ee 

Im =: AGA SEA TEMPERATURE MEA SUREMENT Wa THOULT ATMOS eR IC Ac 2 =NU2 GION 


ime : DIFFERENCE BETWEEN ACTUAL AND AGA SEA TEMPERATURES 


* ALL TEMPERATURES ARE IN °C 
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eo we 27 
SEA SURFACE TEMPERATURE AND EFFECTIVE AT “G" 


DATE eae Soe NOV 27258 = 18: 02 


Wijmesee2e, PeREG-ION = 1.5 m/sec (B-1) 22¢° 

AIR TEMPERATURE 1§.$9°C 

AMBIENT TEMPERATURE Ses. 6 

SEA SUREACE TEMPERATURE ies ; 
Panees DISTANCE : 750m 

Bramintes Guclik ViSTS-o22. 5 NM 

RUMOSPHERIC TRANSMITTANCE : 0. 8842 


BNi@ee eos INCIDENG= AND SUREACE CONDITIONS 


ge wy ese so’ || = og 





i | Nl . 
MEELAT Bet PlEnat | B-4 | |FLAT Sal nnagel 5-4. Tuga Psq< 








| ers | , | | | ! 
R |i¢.804/0.2110/|0.513!0. 138|/0. 29210. 166!/0.180/0.252|lo.:20/¢. 120 | 





|o. €20|0. 849||0. 8680/0. 880 | 


1-R lio. 2196/0. e90/ |0. 48710. 862/ 0. 70810. 834 





15.58/14. 25|/15.06/14.31/|14. 63/14. 74 lag. 40/24. 39} |24.27/24. 27 


eee eetee |e. 25| |=. 06|-0. 31 


| 
| 





| { | | | 
Pecucmenicl Ween 201-0. 49) |=cs 27120127 | 











R Shoe See ANC OF THE SEA SUREACE 

oie Sele eeclm SUREATE 

B-<¢ San SURZACE ROUGHENED BY SEAUFORT 4 WIND 

as CALCULATED APPARENT SEA SURFACE TEMPERATURE 

me STS2SS32NC= BETWHSN ACTUss AND SALCULAcS> SEA TEMPERATURES 

im Rowe Seinen lR= MonSuneMoNT WecHnoUs ATMOS===RIC ATTENUATION 
and Piece neniee SEtna=N ACTUAL AND AGA SEA TEMPERATURES 
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TAB IE 2s 
SEA SURFACE TEMPE RARER D Eee Gree vag tt 


DATS / [iM : 18 NOV 0&:34 - 08:36 
WIND SPEED / DIRECTION : 2.0 m/sec (B-2) 325° 


= = « © 
AIK TEMe snare s See. eS 


AMSIENT TEMPERATURE : 35 ¢ F 
— . —=_= -_= —_— = -_ = ° 
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Pelee iy Ke 
IN-BAND FLUX CALCULATION USING TI-59 PROGRAM 


1. "GENER E 

The program for the T!I-59 calculator used to calculate the in-band flux makes 
use of the Srmpson’s approximation program included in the master library module. It 
Was Originally written for the HP-67 calculator by Robert Pitlak ( Appolo Lasers, 6357 
Arizona Circle, Los Angeles CA., 90045 ). and has been translated for the T[-59. 

To use this program the following values have to be inserted: 

The initial and final wavelengths of the spectral band. 
b. The temperature of the object (in °K or °F or °C). . 


c, ine te number of intervals to be used by the Simpson’s approximation of the 
integral. 


ao : ? 2 
The calculated flux 1s given in Watts.cm* or photons’cm*“sec 


a FORMULAE AND CONSTANTS 
As we discussed in Chapter VI, Section E, Planck’s blackbody law can be written 


in the form 





ho 
O =F. photons'cnisec 
¥ 
ho 
or ies J Cy (q')) dx Watts’cm? 

hy 

l C5 

Where g = ——— x= 
ne (e* 1) jell 


Co = 2mc = 1.883651536 X 10° ( jtm>yem~sec ) 
2mhe* = 37418.42875 ( Watts jim? 'cm? ) 
ch;k = 14388.32334 ( um “K ) 


iat) 
> — 
ll 


The integrations are done numerically using Simpson’s approximation. 
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4. 


PROCEDURE 


How to use the program 


Load the magnetic card 

ness i 

Enter +, initial wavelength in am, press R S 
Enter A> final wavelength in wm. press R S 
Enter temperature (T). press 

Ub tee 

2.8 ifin F 


jc ifin C 


Enter even number of integration intervals 
bess 

= ? 
oe) toes ine answer i 7 atts,cm- 


; a 
(eee eno Cet the ansimemin photons/cm-sec 


PROGRAM LISTING 


In the following Table 24 we give the listing 


calculator. 


of the program used in the TI-39 


TABLE 24 
T]-59 PROGRAM LISTING 


tis M2 2 Lie ne clea 
025” Gia ieee li ee ec 
nS “Gees [ees Bee 175 OS OR 
a) ee liv 4 4 Via) Sami 
as eee Go Lis Soe li - U2 Se 
Qae62 [sa OS? Ue eer Sr ip l, 2? Geo 
De ra bee Hs C2) res Leo) ieee Po) eee li 
ite ee Met os CEL Lil Sse Po eee 
YOS 2 wes leew oe [22> Sere Poe 3S eee 
ws $2 370 i535 “220200 Vee Oise L533 33 AD 
igs 2 ee N54 AN ANZ l2=.yo- see SL es 
i656 22 PRT ee el 22 ee. tou LSE 
OO, Sf RS Une 0 es i Seo eee ned ae 
mn TA LEL eee) ee Loe lie ls, OSS OB DY 
m2 tf Ne> Wo Q bes <)> 12% $2 STO 
O10 Sosa se et) ee Lo? a=) Loan on 1S 
Oi 1° epee Deo Wi ee Lou) Sos. Pa 2 SF 
iS Ses el) Ses ae Lol sex Toi 3s aoy 
Q15 RBS »5 Qv2 O4+ 4 (le 4 Seer Las. Besaaec 
a+ Fee Dt 2 los se 92 eee 
he ani a) Aes. DOME jae ee Caria oy) L3=- NZ ge 
Lilley tel lee) t,o (Bese yoo) SH oo. > aS 
Ole) 808 N75 DA 35 {o5 6fF 8 LBL L325 63 be 
ales ae) yale Nv. o> OF for — 45 2p 122 "OPN oi 
NL? O4+ Os no BE H {3 2S°RDY Ls Soe 
N20, Vereeric Bn ess REL bs> <2 S70 [23.22 an— 
Neo Se Oe et ee ee) fa OS" Bos ony Sec 
Oss 5 2265r N31 35 + L441 99 POT EGiae Uo nas 
Meo ei Ws Qe Cie 2ce [325 22 AY 202 <2" m= 
NZ+ 9s2° OP 25 ny la oe C AAS asc 
es A Sear ia et is <3 let 45 RIEL <Q4 03 OS 
Jes een: = Bie epee Ls epee 23. 9) eet 
Bye fis eee N35 e210 35 TS - 25 = 2220, 
25 elses N37 OO YO ls7r 42 PIL 207 Gs 3g 
022 @e2e > i eee $2 OL D1 Sis Vile ee 
o20 oo oD Ue eae bee 22) Sees 
1S 2 2S) eee lo. SS rie) leat 
B.S 7 ee Net (22570 bol, Waser Jey =5 
Oe Meee ms> OF OF Loc Sails 2 ee 
2s) Sie 0s ae Loe. 205 eis Of 5 
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